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INTRODUCTION 


The  tectonic  history  of  the  Arctic  Ocean  Basin,  and  in  particular 
the  Canada  Basin  and  the  adjoining  Beaufort  Sea  shelf,  today  remains 
unclear  even  after  decades  of  geological  and  geophysical  investigations. 
Reconstruction  of  the  tectonic  history  of  the  margin  requires  information 
about  the  deep  crustal  structure  of  northern  Alaska  at  the  Beaufort  Sea. 
In  1975  and  1976  a substantial  amount  of  refraction  data  was  successfully 
collected  on  the  Alaskan  shelf  from  an  icebreaker.  This  technical  report 
contains  a detailed  description  of  the  methods  and  techniques  used  in 
the  data  collection  and  analysis.  Interpretation  of  the  data  is  the 
subject  of  articles  to  be  submitted  to  professional  journals. 
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DATA  ACQUISITION 

Location 

Personnel  from  Oregon  State  University  and  the  University  of 
Connecticut  conducted  marine  refraction  studies  in  the  eastern  Chuckchi 
Sea  and  the  western  Beaufort  Sea  to  obtain  structural  and  velocity 
information  on  the  continental  margin.  In  contrast  to  most  previous 
refraction  studies  in  the  region  which  were  made  from  stations  located 
on  the  ice,  these  data  were  obtained  using  standard  marine  seismic 
techniques  during  the  Arctic  Summer.  Nineteen  profiles  were  obtained 
between  Wainwright  Inlet  and  Prudhoe  Bay  in  August  1975  from  the  USCG 
icebreaker  Glacier  and  in  August  1976  from  the  USCG  icebreaker  Burton 
Island  with  helicopter  support.  The  profiles,  ranging  in  length  from 
13  to  75  km,  were  roughly  parallel  to  the  coastline  in  about  20  m of 
water.  Heavy  ice  cover  forced  many  minor  course  changes  and  is  the 
reason  for  the  non-linear  direction  of  individual  profiles  shown  in 
Figure  1 . 

Techniques 

A zone  of  ice-free  or  semi  ice-tree  water  exists  along  the  northern 
continental  shelf  area  out  to  a depth  of  approximately  2000  meters  during 
one  month  of  the  year,  and  existence  of  this  zone  permitted  operation  of 
standard  marine  refraction  techniques  with  sonobuoys  and  explosive  charges. 

The  lines  were  shot  in  the  standard  marine  fashion  by  dropping 
charges  in  the  water  from  the  fantail  of  the  icebreaker,  using  powder 
fuses  and  tilt  table,  while  the  ship  was  underway.  Rotor-mounted  yagi 
antennas  received  the  sonobuoy  signals  which  were  recorded  on  a 4-channel 
tape  recorder.  Two  of  the  four  channels  were  sonobuoy  signals,  one  was 
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a clock  signal  and  one  was  a shot-break  signal  from  a streamer  3 m long 
which  was  towed  close  behind  the  ship.  When  radio  contact  with  the 
sonobuoys  was  lost,  at  a distance  of  about  30  km  because  of  earth 
curvature,  a portable  recording  unit  consisting  of  a radio  receiver, 
an  amplifier,  and  strip-chart  and  tape  recorders  installed  in  a 
helicopter,  monitored  the  sonobuoy.  The  helicopter  remained  within 
radio  reception  range  of  the  sonobuoy  during  the  shooting. 

Expendable  naval  sonobuoys  of  the  type  AN/SSQ  41A  were  modified 
for  extended  time  operation  by  addition  of  dry-cell  batteries  which 
worked  well  in  spite  of  the  cold  water.  Explosive  charges  of  less  than 
10  pounds  were  made  up  of  nitro-carbonitrate  (Nitromon)  in  one  pound 
metal  cans.  Explosive  charges  between  30  and  660  pounds  were  made  up  of 
Tovex  in  30  pound  plastic  bags.  Extra  boosters  (Dupont  HDP-1)  were 
required  for  large  charges  or  for  long  burn  times  in  order  to  insure 
complete  detonation  in  the  cold  water. 

Shots  were  detonated  every  three  minutes  at  a ship  speed  of  10  kt 
which  was  slightly  variable  because  of  ice  conditions.  This  resulted 
in  a shot  spacing  of  about  0.7  km.  During  the  helicopter  operations, 
snots  were  detonated  at  intervals  of  from  6 to  15  minutes  resulting  in 
a shot  spacing  of  1.4  to  4.3  km. 

Sonobuoy  deployments  consisted  of  two  sonobuoys  (SI,  S2)  dropped 
at  the  beginning  of  each  line  and  an  intermediate  sonobuoy  (S3)  deployed 
in  the  middle  of  long  lines.  This  resulted  in  a special  line  numbering 
with  letter  A for  the  first  sonobuoys  and  letter  B for  the  intermediate 
one  (for  example.  Lines  22A-23  and  22B-23) 
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The  satellite  navigation  equipment  was  inoperative  during  the  second 
half  of  the  experiment,  therefore  most  of  the  navigation  was  by  radar  fixes 
to  land  points  at  15  minute  intervals.  The  ice  coverage  ranged  from  0 to 
8 octas  during  the  course  of  the  experiment  and  required  frequent  course 
changes  and  caution  on  the  part  of  the  shooter  not  to  drop  the  charge  onto 
pieces  of  floating  ice. 

Sea-surface  currents  were  not  negligible  and  affected  sonobuoy  drift. 
The  direction  and  magnitude  of  the  drift  was  estimated  by  combining  water 
wave  travel  time  and  navigation.  The  in-line  component  of  drift  ranged 
from  0 to  1.65  m/s  and  was  quite  variable  from  line  to  line. 
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DATA  REDUCTION 

A combination  of  manual  and  computer-aided  manipulations  shown  in  the 
flow  chart  of  Figure  2 transform  the  raw  data  to  record  sections  and 
velocity-depth  profiles. 

Arrival  times  for  ground  and  water  waves  picked  on  each  seismogram 
and  other  information  such  as  bathymetry,  ship  velocity  and  streamer 
length,  constitute  an  input  data  file  for  computer  program  TIMCORM. 

This  program  computes  corrected  ground  and  water  wave  arrival  times  to 
a datum,  making  corrections  for  the  shot  instant  due  to  separation  of 
shot  and  streamer  and  surface  and  bottom  corrections  at  receiver  and  shot. 

These  results  form  an  input  file  for  the  computer  program  REFPLTT 
which  produces  a travel  time  plot  with  the  corrected  arrival  times  of 
the  ground  waves  at  the  distances  calculated  from  corrected  water  wave 
arrival  times.  A first  attempt  is  then  made  to  interpret  seismic 
velocities  and  the  data  examined  for  possible  errors. 

From  the  corrected  output  file  obtained  from  program  TIMCORM,  the 
computer  program  REDPL0T2  produces  a reduced  time  plot  where  arrival 
times  reduced  with  a velocity  of  5.00  km/s  are  plotted  versus  distance 

A profile  which  includes  some  helicopter  refraction  data  requires 
an  additional  routine.  A data  file  is  prepared  from  the  ship  navigation 
as  an  input  to  computer  program  UTMSTEVE  and  its  two  subroutine  programs 
UTMGRIDB  and  UTM.  These  programs  plot  the  navigation  data  on  a Universal 
Transverse  Mercator  projection.  Shot  positions  are  added  manually  by 
interpolating  time  and  arcs  are  drawn  using  a compass  at  a radius 
corresponding  to  the  water  wave  travel  time.  Sonobuoy  drift  due  to  wind 
and  currents  appears  as  a non-coincidence  of  the  arcs.  A drift  rate 
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Figure  2.  Flow  chart  of  data  analysis 
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parallel  to  the  trackline  is  assumed  to  apply  to  more  distant  shots  where 
water  wave  arrival  times  may  not  be  detected  by  the  sonobuoy.  The  shot 
distance  is  then  the  distance  between  the  shot  location  and  the  drift 
position  of  the  sonobuoy  at  the  time  of  the  shot. 

The  drift  distances  between  the  shot  and  the  sonobuoy  and  the  ground 
wave  arrival  times  from  the  helicopter  refraction  data  constitute  an  input 
file  for  program  BEAUHELI.  Computer  program  BEAUHELI  computes  time 
corrections  between  the  ship  and  the  helicopter  clocks,  the  corrected 
distance  between  the  shot  and  the  sonobuoy,  and  the  corrected  travel  time 
for  the  ground  arrival.  The  strip  chart  records  played  back  from  the 
helicopter  recorder  are  then  aligned  to  the  corrected  ground  arrival  times 
as  plotted  on  the  reduced  time  plot  and  traced  to  produce  the  record  section. 

Apparent  velocities  and  intercept  times  are  obtained  from  this  record 
section  and  put  into  the  computer  program  MULTLAY.  The  computer  program 
MULTLAY  computes  a model  composed  of  N dipping  plane  layers  from  reversed 
seismic  refraction  profile  data  using  the  formulation  of  Adachi  (1954). 

The  program  requires  two  sets  of  apparent  velocities  and  intercept  times 
from  single-ended  refraction  data.  Complications  such  as  non-reciprocal 
data  or  early  intercept  times  require  a manual  ray  tracing  method  based 
on  the  formulation  of  Adachi  (1954). 

On  each  seismic  model  obtained  from  MULTLAY  the  portion  of  an  interface 
drafted  with  heavier  lines  on  Figure  8 through  Figure  26  represents  only 
that  portion  of  the  layer  which  gives  rise  to  observable  seismograms.  The 
offset  distances  of  the  rays  are  obtained  from  the  program  WAREFRA.  The 
computer  program  WAREFRA  is  based  also  on  the  formulation  of  Adachi  (1954) 
and  requires  the  seismic  model  obtained  from  MULTLAY  as  an  input  data  file. 


For  the  combined  velocity-depth  section  (Figure  27),  each  profile 
was  projected  onto  a straight  line  drawn  through  the  tracklines. 


A. 


DATA  INTERPRETATION 


Because  of  the  navigation  difficulties,  only  five  lines  were  reversed. 
These  lines  are  Lines  1-2  and  6-7,  Lines  3-4  and  4-5,  Lines  10-11  and  12-13, 
Lines  22A-23  and  24B-25,  and  Lines  24A-25  and  22B-23  where  the  numbers 
indicate  the  end  points  of  the  lines  (see  for  example  Figures  10  and  11). 

The  nine  remaining  lines  (Lines  GL-1,  GL-2,  8A-9A,  8B-9B,  16-17,  18-19, 

20-21,  26A-27A,  26B-27B)  were  interpreted  as  single-ended  lines  (see  for 
example  Figure  19).  For  complicated  lines  showing  non-reciprocal  data 
or  early  intercept  times,  a manual  ray  tracing  method  was  used  (see  Figures 
16  and  17). 

Figures  8 to  26  show  the  nineteen  profile  plotted  as  reduced  record 
sections  where  the  reducing  velocity  is  5.0  km/s.  The  distance  is  plotted 
as  water  wave  travel  time  seconds  where  the  water  velocity  was  estimated 
to  be  1.44  km/s.  All  the  velocities  are  in  km/s.  The  seismic  models 
determined  from  interpreted  reduced  sections  are  composed  of  plane  horizontal 
or  plane  dipping  layers.  These  are  shown  above  the  record  sections  with  a 
vertical  exaggeration  of  3:1.  Heavy  lines  on  the  refractors  indicate  the 
interfaces  responsible  for  observed  arrivals  and  contain  the  horizontal 
offsets  for  upgoing  rays.  Only  the  reversed  profiles  resulted  in  true 
velocities,  and  on  all  the  single-ended  profiles  apparent  velocities  are 
assumed  to  be  true  velocities.  A few  estimated  velocities  are  shown  in 
parentheses. 

The  composite  section  shown  in  Figure  27  summarizes  the  subsurface 
velocity  and  depth  information  calculated  from  the  refraction  data.  The 
heavy  lines  on  the  interfaces,  corresponding  to  the  observed  arrivals  as 
shown  in  Figures  16  through  26,  were  projected  onto  a composite  section 


which  passes  through  the  profiles.  The  small  arrows  along  each  interface 
correspond  to  the  appropriate  refraction  lines  shown  at  the  top  of  the 
figure.  The  refraction  layers  from  adjacent  lines  have  been  correlated 
on  the  basis  of  velocity  as  indicated  in  Figure  27  by  light  lines. 
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2 . C 

3 . C 


3 . C 


UTN$  T£v£ 

COMPILER  DOUBLE  PRECISION 

PROGRAM  UTMMAP 

7 h I $ PROGRAM  CALLS  UT  MGR  19  70  *arE  a N UNIVERSAL  7 p * .*»  S V £ R 3 E 
N E *C  A T 0 S GRID  IT  T i i E N R £ h i j mH  AES  0 MAG  DATA  FILE  * N D P L - T S 
EITHER  TSACKLINES  OR  AMNGThTED  DATA  POINTS  OH  THE  GRID 
COMMON  /UTM2/  NLAT.  SLAT.  EJ.3MG,  ALONG 
REAL  LAT.  LONG.  NlAT,  LX,  If.  NX.  NY,  LATH,  LON  GY 


IB. 

PI  ■ 3 1419924934 

Hi 

P 1 92  • PI  ' 2 % 

1 2 . 

P 1 9 4 • PI  / 4 i 

13, 

TWOPI  • PI  ♦ PI 

14  . 

R TO D • 19B  / PI 

19, 

9 TO  R • PI  / 18B 

14. 

tSl9E  - l 

17,  C 

GO 

CET  THE  CBIB  PABANETEBS  AHt  DBA*  THE  CBIO 

IB, 

CALL  UTMCRID 

IB. 

TYPE  •ANNOTATE  THE  DATA  PT$' 

28 , 

READ  (11.  1 8! > IANS 

21  . 

IBI 

FORMAT  <A2) 

22.  C 

SET 

UP  THE  INPUT  FILE 

23. 

9 

TYPE  •FILENAME- 

24. 

READ  (11.  1 92)  FNAME 

' I ) 

29. 

1B2 

FORMAT  (SIB) 

24  . 

CALL  FQPEN  (l.  FNAME) 

27.  C 

GET 

A DATA  POINT 

2«. 

IB 

READ  (1,  IBB.  END-28) 

I VAR  I,  LAT.  LATM,  L 

29 . 

IBB 

FORMAT  ( 7X.  2 A 2 . IX,  F3 

9.  IX,  F9  9.  IX.  F 4 3.  IX 

1B.C 

CHECK  LPT  AMI  LONG  FOB  HITHIN  THE  G» I 9 

Jt. 

LAT  • LRT  ♦ SIGN  (LATM/48  . LAT) 

32  # 

LONG  • LONG  ♦ SIGN  CLONCM/48  . LONG) 

33. 

IF  (LAT  EO  t •>  GO 

TO  IB 

34. 

IF  (LONG  LT  B B > LONG  » 343  * LONG 

39, 

IF  (LAT  LT  SLAT  OR 

LAT  GT  NLAT)  GO 

34. 

If  (LONG  LT  WLONG 

OR  LONG  GT  ELONG) 

37.  C 

CONVERT  LAT  AND  LONG  TO 

UtM  X AND  Y 

3B . 

CALL  OTM  (LAT.  LONG. 

X.  Y) 

19  . C 

IF 

NO  ANMOATAIGN  DON'T  NEED  TmE  ANGLES  ETC 

*8. 

IF  (IANS  EO  2NN0)  CO  TO  14 

*i  . C 

FIGURE  OUT  tfMCRC  AH 9 WHAT  ANGLE  TO  ANNOTATE 

42  . 

ISI9E  ■ ISIDE  • (-1) 

43. 

ynly  * r - ly 

44  . 

THETA  • ATAN  (YMLY  / 

<X-LX>> 

49. 

LX  • X 

44  . 

LT  • Y 

« • 

IF  <A8$  (THETA)  lT 

•134)  GO  TO  12 

48  . 

ANGLE  • B 

4*. 

NX  ■ X 

9B 

IF  ( ISIDE  LT  9)  NX 

» NX  - 9 3 

91  . 

NY  « Y 

92 . 

GO  T0  13 

93. 

12 

ANCLE  • 2 78 

94  . 

NX  • K 

93. 

NY  • Y 

96 

IF  f ISIDE  LT  3 > NY 

* NY  ♦ 8 3 

97.  C 

LPOEl  THE  DATA  POINT 

31 

19 

CALL  SYMBOL  (NX.  NY, 

B 87.  IVARI.  ANGLE. 

3f  ; 

CALL  PLOT  (X.  Y.  3) 

4B  . C 

PLOT  THE  IATA  POINT 

€1  . 

14 

CALL  PLOT  (X.  Y.  2) 

42. 

CALL  MARKER  ( 1 > 

€3. 

GO  TO  IB 

44.  C 

HEBE  ON  EOF 

49. 

2B 

CONTINUE 

44. 

CALL  FCLOS  (1. 

47  . C 

ANOTHEB  FILE  ONTHIS  SB  ID 

43. 

TYPE  'ANOTHER  FILE" 

49. 

READ  (11.  1 9 l > IFILE 

LONGM 


78. 


If  £iriL£ 
END 


EO  2HYE)  GO  TO  5 


I 
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KTNdill 

conduct  BOUtLC  FRCCISION 
SUMOUTINC  U 7 NCR  X | 

THIS  SUBROUTING  MAMS  A UTN  NCRCATOR  GRIB  ON  THC  FLOTTCR  17 
RUCRR1CS  7NC  U SC  R FOR  TNf  NAP  PARANC  7 |RS  AN  | HAVES  TNCN  IN 
CONNON  FOR  7 NC  CALL  INC  PROCRAN 

S0U7NCRN  LA  7 I 7 U II S *N|  MCS  T CRN  L ONC I T UBCS  ARC  NfC A 7 1 V| 

U7NNAF  MILL  MORI  RCROtf  R07H  7NC  IA7  CL  I N|  ANl  CRCCNMICN 
NCR1B1 AN 

RIAL  NINLA7.  NAXLAT.  NINLON*.  NAXLONG 

CONNOR  /U7N|/  SCALC.  XO.  TO.  9T0R.  R70I.  A.  R.  C.  LA.  LALA. 

I L ISO • ONLfSO.  CNCRIB.  10 

RIAL  LA.  LO.  Lf  SO.  LALA.  NTOI.  10 
RIAL  NINLATN.  NAXLA7R.  NINLONCN.  NAXLONCN 

UAL  LAf.  LONG.  NL  AT.  NIN 

CONNON  /UTN 2/  NLA  7 . SLAT.  (LONG.  ML  ONC 
S ININS  ION  AiCI  <40) 

COBCS  FOR  TNC  FL  0 77|R  SU9R  00  7 INCS 

1NTCCCR  FCNUR.  FCNSOMN.  PLOTCNB.  FIN.  ARCD 
FCNSOMN  • 2 
FCNUF  • 3 
FLOTCNS  • -FCNUF 
C0NSTANS7  FOR  CAL CUL A 7 IONS 
•TOR  • 3 141392(931  / ISO 
N 70  X • 39  IF 

PARANCTCRS  for  a utngris 
10  • 0 9990 
A • I •091091 
t » t.  tismt 
C • 0 0000100 
LA  • S3F02BS  A 
LO  • 09(911  0 
LALA  • LA  • LA 

LCSi  • <LALA  - LO  • LO)  / LALA 
0NLC90  *10-  LCSO 

XO  • 0 0 
TO  • 0 0 

iS  MR! TC  (10. 100) 

LONGS  MIST  ARC  NCCATIVC  ANB  LATS  SOU  7 N ARC  NIC A 7 1 Vl 
ISO  FORNAT  < SGIVC  7N|  FOLLOMX NG  ' . '. 

1 ' NL  A 7 SLAT  MLONG  CLONC  SCALC  Gill''. 

2 BNSNSNINX  N ) 

RCAS  <11.  1 0 1 > NAX'.AT.  NAXLA7N.  N 1 NLA  7.  NINLATN, 

1 NINLONG.  NINLONCN.  NAXLONG.  NAXLONCN. 

2 SCALC.  NIN 

1 V I FORNAT  < F 3 A.IX.F2  R.1X.F3  B.1X.F2  O.IX  Ff  B.iX.F?  i.lX.rA  A 

t ix.  PS. o. «x. re  o. ix. fa  e: 

NAKC  ICCINAL  DCCRCCS  OUT  OF  T NC  SCGRCCS  ANB  NINU7CS 
SLAT  • N X NL AT  ♦ SIGN  (NINLATN.  NXNLAT)  / SB 
NLA  7 • NAXLAT  ♦ SIGN  (NAXLA7N,  NAXLA  T ) / SB 
ML  0 NC  • NINLONG  ♦ SIGN  (NINLONCN.  NINLONG)  ' SB 
CLONC  • NAXLONG  ♦ SIGN  (NAXLONCN.  NAXLONG)  / SB 
NIN  • NIN  ' SO 

IF  ( NL AT  LC  SLAT)  CO  TO  IS 
N All  LONGS  0 TO  3S0  TO  TN|  CAST 

IF  (MLONG  LT  • ) MLONG  • 3SB  * MLONG 
IF  (CLONC  LT  • > CLONG  • ISO  ♦ CLONG 
T All  CAR!  OF  CROSSING  CRCCNMICN  - NAKC  LONGS  GO  TO  ?2B  DIGS 
IF  (MLONG  LT  CLONG)  CO  70  12 


• J 

J 

f 


■ 
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(LONG  • C LONG  ♦ 36* 

6*.  12  MRI TE  ( It.  I \?  ) 

<i.  112  format  <-  garret  read*  • 

62 . RE  AO  (11.  114)  IANS 

63  114  FORMAT  <A2) 

64.  IF  (JAMS  CO  2NYC)  GO  TO  19 

<9 . GO  TO  1« 

66.  16  TYFE  * ML  A T . SLAT  ERROR* 

62.  GO  TO  1R 

69  12  TY»C  'LONGS  NOT  Ml T M I N • OR  - 3 DEGREES  OF  CENTRAL  MERIDIAN' 

69.  GO  TO  l« 

**.  19  CALL  F OF  C N <6.  'IRLT*) 

21.  C INlTiaLUg  the  flotter 

22.  CALL  INITIAL  (6.  1 ft.  -•  9.  22  > 

73  . CALL  FlOT  ( l . 1 . FLOTEND 

74  C INITIALIZE  FEN  POSITIONS 

'9.  SCALE  • l t / SCALE 

76. C LETS  LABEL  IT  AGAIN  mITN  INFORMATION  CONCERNING  'HE  NAP 
* 7 . C AREA.  DATA.  ETC 
. MR  l TE  (It.  1 1 2 > 

*9  1 1 2 FORMAT  C'tcivc  TNC  PLOT  LAtEL  --  UP  TO  It  CHARS  > 

St.  READ  ( 11 . 1 1 3 ) At C D 

91  . It3  FORMAT  < 4 t A 2 > 

1 2 . CALL  SYNtOL  (t  -1  . 2 1 . At C 9 t tt> 

9 I . TfNF  - ML  ON  G 

64.  IF  * ML  ON  G GC  1ft  > TtNF  • MLONG  - 36« 

89.  I ZONE  • (TEMF  ♦ 116  > ' 4 

86  CMERID  • NOD  ((IZONE-l)  • 6 ♦ 183  369  ' 

67.  C NOTE  CNERID  9 TO  369  SCGtCES  TO  Tm£  EAST 

99.  IF  * DAf S (CMERID  - ML  CNG  > CT  3 'CO  r0  (7 

39.  IF  iDARS  (CMERID  - ELONG)  GT  3 > CO  TO  17 


«t . C 

NOM  Ml  TNI N ♦ 

OR  - 

3 DEGREES  OF 

CNCRI D 

31 

IF 

< MLONG 

GC 

CMERID) 

CALL 

utn  (Slat. 

•LONG. 

92 

IF 

v CL  3NG 

LE 

CMERID) 

C ALL 

UTN  • SLAT 

ELONC. 

9 3 . 

IF 

( MLONG 

LT 

CMERID 

AND 

ELONG  CT 

C NCR  I D > 

34  . 1 CALL  UTN  (SLAT.  CNERtO.  <.  Y) 

?•  . vt  ■ Y 

36. C NOM  HAVE  N I N Y VALUE  ( Y0 i FOR  THE  NAP 

37  IF  (MLONG  GT  CMERID'  CALL  UTN  CNLAT.  MLONG.  y f ) 

?9.  IF  t ML  CNG  LE  CMERID*  CALL  UTN  SLAT  VLDNC  X v> 

99  Xt  • X 

1 9t . C NOM  NAVE  MIN  X VALUE  v Xt  ' FOR  THE  MAP 

191. C <49.  Yt)  are  the  SOUTH  m£$t  CORNER  OF  The  map  Such  ’**<»•  *-£  ; 9 : “ 

142. C DvCC  NOT  GO  E l THE  P SOUTH  OR  MEST  Or  THIS  POINT 

1*3. C GUESS  that  me  CAN  START  MAKING  THE  GRID 

104. c START  the  grid  DRAMING  the  lats  EAST  and  MEST 

1*9.  LONG  • MLONG 

106.  LAT  • SLAT 

102  DLONG  • (ELONG  * MLONG)  ' !•  *t  1 

If*.  CALL  UTN  (LAT.  MLONG.  X.  V' 

199.  CALL  NljMtER  <X-  79.  V.  14  N I N L A Y 9 - 1 ) 

lit.  CALL  NUMBER  (X-  39.  Y.  14.  MINLATM.  9 . -1) 

III.  CALL  FlOT  < X.  Y.  FENUP) 

1 12  . FEN  • PE  N DO  MM 

III.  CO  TO  29 

114. c DRAMING  A LAT  FROM  MEST  TO  EAST 
119.  29  FEN  • PCNUP 

116.  LAT  • LAT  ♦ MIN 
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1 17  ■ 

DLONG  • -DLONG 

119 . 

LONG  • ML  ON G 

1 1 9 t 

IF  (LAT  GT  MLAT)  GO  *0 

I It. 

23 

IF  ( LONG  GT  EL  OH C ) 00  TO  39 

Ill  » 

CALL  U TH  < L AT . LONG.  X.  Y> 

122. 

CALL  PLOT  (X.  Y,  PEN) 

123. 

PEN  • PEMDOVN 

1 24  . 

LONG  • LONG  * DLONG 

lit. 

GO  TO  29 

129.  C 

BRAVING  A LAT  FRO*  THE  EAST  TO  VEST 

127. 

39 

PEN  ■ PENUP 

129. 

LAT  • LAT  ♦ MIN 

129* 

IF  (LAT  GT  NLA T ) GO  TO  4$ 

1 39  i 

LONG  • ELONG 

131  . 

DLONG  • -9L0MC 

1 72) 

CALL  Um  (LAT.  LONG.  X.  Y > 

1 33  . 

CALL  PLOT  (X.  Y.  PEN) 

1 34  . 

PEN  • PEN  DO VM 

t 39  . 

33 

LONG  • LONG  4 DLONG 

1 39. 

IF  (LONG  LT  NLONG)  GO  TO  29 

1 37  i 

CALL  U TH  (LAT.  LONG.  X.  Y) 

1 39. 

CALL  PLOT  (X.  Y.  PEN) 

1 29  ; 

GO  ?0  33 

1 49. 

49 

CONTINUE 

1 4!  . C 

tout  VfTH  EAST  AMO  v£ S T PORTION  Of 

THE  GRIB 

M2.  C 

MOM  LETS  BO  THE  NORTH  AMI  SOUTH  LINES 

1 43  . 

LONG  • MLONG 

1 44  . 

LAT  • SLAT 

1 43  . 

DLAT  ■ (NLAT  - SLAT)  / 19  991 

1 4C  . 

CALL  UTN  (LAT.  LONG.  X.  Y> 

14  7 i 

CALL  NUH9ER  (X-  49.  Y - 23.  14. 

H I ALONG. 

1 49  . 

CALL  NUH9ER  (X4  it,  r-  23.  14. 

H I MLONGN. 

1 49. 

CALL  PLOT  (X.  Y.  PENUP) 

139. 

PEN  • PENBOMN 

191  . 

GO  TO  29 

1 72.  C 

CRAVING  A LONG  FROM  SOUTH  TO  MORTH 

133. 

43 

PEN  • PENUP 

134. 

LONG  ■ LONG  ♦ H1N 

! 33) 

DLAT  • -OLAT 

1 76  . 

LAT  • SLAT 

1 27. 

IF  (LONG  GT  ELONG)  GO  TO  79 

129. 

39 

IF  (LAT  GT  NLAT)  GO  T 0 33 

’.29  ; 

CALL  U TH  (LAT.  LONG.  X.  Y) 

1 S9. 

CALL  PLOT  (X.  Y.  PEN) 

191  . 

PEN  ■ PENDOWN 

» 92  i 

LAT  ■ LnT  4 DLAT 

1 93  . 

GO  TO  39 

194  . C 

CRAVING  A LONG  FROM  NORTH  T0  SOUTH 

1 93  t 

33 

FEN  • PENUP 

1 99  . 

LONG  • LONG  4 HIM 

197  . 

IF  (LONG  GT  ELONG)  GO  *0  70 

1 99  . 

LAT  • NLAT 

1 99  . 

OLAT  - -OLAT 

179. 

CALL  U TH  (LAT.  LONG.  X.  *> 

1 71  ; 

CALL  PLOT  (X.  Y.  PEN) 

1 72) 

PEN  • PENOOUH 

1 73. 

99 

LAT  • LAT  ♦ OLAT 

1 74 

IF  (LAT  LT  SLAT)  GO  TO  49 

1 79  . 

CALL  U TH  (LAT.  LONG.  X.  t> 

1 79. 

CALL  PLOT  (X,  Y.  F£N) 

1 77. 

GO  TO  99 

1 79.  C 

DOME  9ZTM  THE  GRID 

i 79. 

79 

CONTINUE 

199. 

RETURN 

191  . 

ENI 

l i 
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uT« 


1 . 

CO.. Il£.  30USLE  .8ECISIOA 

2. 

suo.odtihe  uta  clat.  lons.  x.  t> 

2.  C 

this  suoaoutine  contexts  lat  (.hd,  long  '3laao»a>  To 

X. 

Y 

4,  C 

IN  J TA  C0083INATES  IN  INCHES  E80A  THE 

WOOER  LEFT  CORNER 

OF 

3.  C 

THE  .LOT  AT  SCALE 

5. 

SEAL  L«.  L 8 . LE  SO  < LA-  NU.  LALA. 

HTOI.  LAT 

LOMQ. 

<0 

7 . C 

jTAt  HAS  .OOjECTIO*  .A.ANfTE.S  That  AtE  >AS3E3 

*»0H  UTHGRI3 

9. 

CONNON  'UTAIt  SCALE.  XO.  TO.  ST08 

. HT01.  A 

8.  C. 

LA 

o LAIh, 

9 * 

1 LESS.  OALESO.  CAEBU.  ICO 

19. 

.Ht  • LAT  • 0 T 00 

11  . 

SIN. MI  • SIA  CHI) 

12* 

S.S.  • SIN.N1  • SIN. HI 

13. 

OHO  • LA  • OALESO  / <1  a - Ltsa  • 

$PSP>  •• 

t 3 

14  . 

hu  • la  / satt  ci  a - L£sa  « s.s.> 

13. 

LA  • LA  • OALESa  • (A  • .NT  - 0 • 

SIN  (2  a 

PHI)  / 

2 

♦ 

I C • SIN  < 4 • .HI)  / 4 > 

17  i 

3L  A AID  A • ( LONG  - CAE8I3)  • 3TOO 

19. 

3L»L  • OLAAODA  • OLAASDA 

19. 

OLSLlL  • DL  DL  • DLAAttA 

19. 

SIN. HI  • SIN  (.MI) 

21 . 

COS. MI  • COS  (.MI) 

22. 

TAN. HI  • TAN  (.HI) 

23. 

y • K9  • <LH  ♦ OLOL  • HU  / 2 9 • SlNPHl  • COSFHl) 

COSPM  1 

4 ♦ 

X a xf  • (0LAH99A  • HO  • CQSPXl  ♦ 

9L3L3L  • 

HU  / 9 

8 • 

23. 

1 <HU  ' RmO  - 7A HPH1  • TAHPHl>> 

29.  C 

SCALE  X A H 1 Y TO  IHCHES  OH  THE  PAPER 

27. 

X • X • SCALE  • HTOl  - <« 

29. 

Y a Y • SCALE  • HTOI  - Y9 

29  . 

tCTURN 

29. 

EH9 
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24 

20 
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30 

69 
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71 

94 
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47 

50 
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33 
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71 
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71 
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23 
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43  -134 

39 

99 

31  760011 

1990 

71 

13 

63  -134 

24 

39 
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7 1 
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44  -134 

19 
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00 
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00  -153 

27 

80 

42 , 7608  1 1 

23  39 

7 l 

0 1 

10  -153 

22 

39 

43  7608  1 1 

2343 

7 1 

0 1 . 

10  -153 

18 

30 

44 . 7698  1 2 

0908 

7 l 

0 1 

56  -’33 

13 

00 

43  760812 

09  • 3 

7 1 

0 1 

70  -153 

e: 

49 

46  766812 

0039 

7 1 

00 

60  -132 

58 

99 

4 7 . 7608  1 2 

0943 

7 l 

99 

99  -132 

35 

29 

48  7608  1 2 

0100 

70 

39 

39  -132 

47 

30 

49, 760812 

01  1 3 

70 

59 

29  -132 

40 

60 

30 . 7608  1 2 

0130 

70 

38 

39  -132 

33 

40 

31 . 760612 

9143 

79 

38 

70  -132 

39 

89 

32. 760812 

0290 

70 

59 

39  -132 

43 

49 

33  « 7609  1 2 

0213 

70 

39 

89  - 1 32 

59 

49 

34  . 7698  1 2 

9230 

71 

00 

59  -132 

54 

39 

33. 760812 

0243 

7 1 

91. 

39  -133 

09 

90 

36. 760812 

9390 

71 

91 

39  -133 

92 

29 

57 , 760012 

0313 

7 l 

• 2 

19  -133 

•6 

90 

30. 760012 

0330 

71 

• 2 

19  -133 

12 

70 

39. 768012 

0343 

71 

02 

79  -133 

19 

00 

60,760912 

0400 

71 

02 

99  -133 

24 

00 

61 , 760912 

0749 

71 

99 

99  -133 

92 

40 

62. 760012 

• 999 

71 

90 

49  -132 

53 

30 

63. 760ei2 

9020 

79 

39 

33  -132 

44 

60 

64, 760012 

0900 

70 

38 

63  -132 

39 

33 

63 • 760812 

0922 

79 

38 

•9  -132 

32 

30 

66. 760912 

0949 

79 

36 

93  -132 

23 

99 

67. 760012 

1999 

79 

33 

63  -132 

17 

29 

60. 760912 

1920 

70 

33 

90  -132 

09 

39 

69. 760812 

1040 

79 

54 

29  -132 

99 

70 

70. 768812 

1190 

70 

33 

43  -131 

33 

99 

71 . 760012 

1129 

79 

33 

29  -131 

43 

39 

72. 760912 

1149 

79 

3 1 

19  -131 

33 

99 

?3.7S*«12  12**  7*  ;e  2*  -HI  2«  «• 

74  76*812  1221  7*  4*  5*  -131  13  «< 

73.76**12  124*  7*  47  3*  -131  1*  *• 

76.76*812  13**  7*  47  »•  -131  *3  8* 

77.76*812  132*  7*  46  3*  -138  33  8* 

71.76*812  134*  7*  44  7*  -138  48  88 

73.76*812  14**  7*  44  ••  -13*  4*  88 

8*. 76*812  142*  78  41  42  -138  32.3* 

81.76*812  144*  78  41  32  -138  21  38 

*2.76*812  13««  78  4*  88  -13*  1*  4* 

83.76*812  132*  78  48  18  -138  82  8* 

84.76*812  134*  78  38  3*  -143  32  3* 

83.76*812  168*  78  37  68  -143  43. •• 

*6.76*813  *34*  7*  32  12  -148  33  3* 

*7.76*813  *6«*  7*  32  S3  -148  46  8* 

•* i 76*8 1 3 *628  78  34  66  -148  36  28 

*3.76*813  8648  78  33  78  -143  83  48 

3* i 768813  87*2  78  33  -143  13  3* 

31.76**13  *722  7*  2*  -143  23  t* 

32.76**13  *812  78  33  8*  -143  3*  8* 

33.768813  *82*  78  36.8*  -143  37  38 

34.76**13  *84*  7*  37  18  -143  3*  73 

33  * 76*813  •*«•  78  37.47  - 1 43  36  •• 

36 . 76*813  *32*  7*  37  23  -IS*  87  2* 

37.76**13  *33*  78  36  88  -13*  23  3* 

3*. 76*813  18*8  7*  36  73  -13*  38  8* 

33.76*813  1*23  7*  37  43  -13*  39  ** 

188. 76**13  114*  78  37  33  -13*  38  68 

1*1.76*813  128*  7*  37  43  -138  27  8* 

1*2.76*813  122*  78  3?  89  -138  16  91 

1*3.76*813  1242  7*  37  18  -13*  83  63 

1 84  * 76  88  1 3 138*  7*  37  4 1 - 1 49  3*  I* 

1*3.76*813  132*  7*  37  3*  -149  37  2* 

1*6.76*813  1341  78  37  3*  -149  37  2* 

1*7,76*813  14*8  7*  37  3*  -149  37  it 

I**. 76**13  142*  7*  37  3*  -149  34  8* 

1*9.76*813  1433  78  36  67  -149  46  83 

11**768813  13*8  7*  36  23  -149  36.23 

111.76*813  132*  7*  36  3*  -149  3*  73 

112.76**13  134*  7*  33  38  -149  28  «• 

113.761813  16*8  78  33  38  -149  89  8* 

114.76*817  *38*  71  8*  18  -132  8*  2* 

113*76*817  *321  7*  39  If  -131  34  41 

116. 76*817  *34*  7*  37  9*  -131  47  1* 

117.76*117  *6**  7*  36  6*  -131  3*  2* 

lit*  76**  1 7 *62*  7*  36  1 * - 1 3 1 27  7* 

1 19  * 76**  1 7 *64*  7*  33  3*  - 1 3 1 1 7 7* 

1 2*  * 76 »• 1 7 *7**  7*  33  9*  - 1 3 1 *8  3* 

121.768817  *72*  78  33  9*  -131  *8  3* 

122.76*817  *748  78  32  98  -131  8*  4* 

123.768817  8*8*  78  31  88  -138  3*  7* 

124*  7688  1 7 *82*  78  4 7 28  - 1 38  38  *6 

123*768817  883*  78  43.88  -138  42  88 

126*768817  8988  78  44  88  -138  39  88 

127  76021 7 092;  70  42  38  -138  24  38 

.28.761317  8338  76  39  39  -13*  12.72 

129.76*818  *388  78  4*  33  -13*  C* . 2* 

13*. 768818  832*  78  41  68  -13*  18  7* 

131.76*818  834*  78  42  60  -138  28  28 

132*760618  *48*  78  44  28  -130  29  88 

133.768818  8420  70  43  88  -130  39  20 

134.760818  *44*  7*  47  3*  -138  48  *• 

133*768818  8308  78  47  4*  -138  38.88 

136*768818  8328  70  49  88  -130  32.0* 

137*760818  8340  70  49  20  -131  03  80 

138*760818  0668  78  30  68  -131  13.68 

139*766618  8620  76  32.68  -131  24  38 

148*766818  864*  78  34  88  -131  33  *• 

141*768818  *788  76  34  86  -131  44  68 

142*760610  0728  76  33  68  -131  33  88 

143*768818  *748  78  36  78  -132  88  8* 

144,76*818  108*  71  82  88  -152  39  4* 

143.768*18  182*  71  84  12  -133  83  It 

146.768818  1648  71  83  66  -133  13  28 

147.760818  lift  71  86  23  -133  26  1* 

148.76*818  1128  71  87  8*  -133  33  2* 

149.76*818  1148  71  87  43  -133  36  3* 

138.768818  1288  71  87  68  -133  42  28 

131*768818  1228  71  88  48  -133  32  88 

132*768818  1243  71  89  17  -134  83  68 

133.768818  1388  71  16  06  -134  14  40 

134.768818  1326  71  II  48  -134  26  S3 

133.766818  1342  71  12  28  -134  38  48 
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COMPUTER  PROGRAM  TO  FORM 
ARRIVAL  TIME  CORRECTIONS  (TIMCORM) 
ANO  INPUT  AND  OUTPUT  FILES 
FOR  ALL  LINES 
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program  n*co*f 

01 Mf NSIOh  G(22l .GCf?2l • AC 111 tR!Af 
»EAC<5,51AI01,AI02 

5 F0RNATC2A81 

6 Rf  AO (5. 101  OBPN,  OMBM.OeEM,  TOT  S.ORM, VI. 3VNS.CLM,SL«.S0M 
10  Fot*AT!3F5.f,7F5.2i 

wRITEfT, 1AI AlOt  «AI02 
1%  F09«*AT(?A8) 

WRI TE ( 7,  1%)  TOTS*  VI 

16  CQtWATC 15I,F5.2.5*,F5.3I 
IS  PO**ATt<  *.2*5, /I 

WRTTE!61.t5> AlDttAIO? 
wRITE(61  an 

17  FORMAT!#  OBRN  O0B»  O0CM  TCTS  OR*  Vt  SV»S  #. 
I#CLN  SIN  SON#) 

WRITE (61  #201  CQRH.OEBN.oeCH.TCTS.ORM.vi  ,SfMS.CLN,SI.M. 
ISON 

20  'OftMAfftf  i ,JF6.0  ,2F6.  2.F6. 3,  AF*  .2) 

WRITEI61 • *** 

25  *ORMAT(#  SHOT  TINE  OISTS  TS  TC  TCOR  Cl  G2#« 
1*  G3  GA  G5  G6*l 
M*0 
ILL*0 

31  ttEACt9,35»ACE'<Af I) ,I*lvll> 

35  FORMAT  ( A5  , 1 1F5  . 0) 

IF(F0F(5II  ILL  *1 

IFfN.EO.ll.ANO.ILL.CQ.il  GO  TO  A0 
IFfUL.EO.il  CALL  EMIT 
IF  UCE.EQ.5MS5599I  GO  TO  % 

IrCACC.EQ.5M  1 GO  TO  135 
I*TFE«0 

AC  IFf  ITTFE.EO.  O.ANO.N.EO.  01  GO  TO  55 
0*0»TC 
0ISTS«0»TS 
00  A3  T*1#N9 
A 3 Q f 1 1 *R  f T I ♦ TC 

WRITE! 7. ASISmot, TT TM • CSSM,TS , 0 • (R 1 1 1 • I >1, NRt 
A5  F0RMATCA5. I5.F5.0.15V.6F6.2I 
00  A 6 T*1.N 
G(I)*GCII*TCCR*TC 
A*  GC I Il*Gt I) ATS 

IFfN.CT.UIGO  TO  55 
WRT  TE ( 7, *0 1 (Gf 1 1 . 1*1 , Ml 
«3  FORMAT  15*.  I 1F5.  21 
GO  TO  50 

55  WRITE «7,50l IGfXl • i*  1 • 11) 

WRITE! 7, 501 !G! II •1*12, Ml 
50  KN*N 

IF(N.GT.6IKN*4 

WRITE (61 ,65  I SHOT .ITIM*  01  STS • TS ,T C. TCO B. ( GC I H , I * 1 . KNI 
IF  fKN.EO.NI  GO  TO  *0 
65  FORMAT ( * t ,A5,I5,10Fe.2T 
KN*7 

70  IFfN.6T.KNA6!  GO  TO  75 
<N*N 

GO  TO  SO 

75  KN*KM45 

SO  WRITEI61 ,851 (GC(I1 •I*KM,KN) 

85  FORMAT!*  *,3AY,6F6.2! 

IF ( KN. FQ .N) GO  TO  50 

KN«KMA5 

GO  TO  70 
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90  IFrrLL.EO.il  CELL  E*IT 
TFlH.EO.m  GO  10  94 
N*0 

00  TO  50 

94  N*0 

99  SHOT *ACF  S niH.llll  | 08S«I(JI  t »2*4<J1 
8T*I»5I  t TS*»16I  10**171 
NR*0 

oo  iso  r*t.4 
RlI)*AtI»7> 

rF(»<II.£0.0.l  GO  TO  109 

too  N**r 

1 09  S1*SQ0T(  1.-V1*V1/VJ/VJI  /VI 
S2*S0RT ( 1.-V1*V1/V2/V2I /VI 

s*=eT*sv»s 

ST=- <.0.99*1. 91  5»9T 

rFisT.GT.oesHiiy.ossn 

ifis*.5T.cl«igo  to  no 

TC*SQRf < (9  T-CL91  ••2*1  SV-SOH)  ••El*.  001/Vt 
GO  TO  120 
110  TL9»CL9*SL9 

IF1SX.GT.TLN1  SO  TO  119 
MV*S0H-SV 

rc*tes(wvi*.soi/vi 

GO  TO  120 

115  TC*SORTMSX-TLM)  ••2*C  ST-SOH)  ••21*.  001 /VI 
126  COHTINUE 

ang=asinfivi/vji 

OFFSET*OBS9»T«NF IlNGI*. 001/VI 

DHS*08*I‘*<08E9-0I’891»  10  *TS*TC  -OFFSET I / TOTS 

offset  =OBRIS*TA  NF  (ANGI  •.  001/VI 

090=08  »HMDBEN-0?BNI  "OFFSET /TOTS 

C*1 09B-08RH*09S-08SH) *. 001 

!Ff v2.nE.VJI  GO  TO  125 

OT1  * 1ST  •0091  »S  2»  . 00  1 

OT  ?=C*S2 

GO  TO  13" 

125  SI=SGRT<  l.-V2»V2/V  J/VJ1  7V2 
0Tl=(ST»C*H1  •.  3 0 1 *S  l 
DT?-.C»CSl-SJI 
1J0  CONTINUE 

TC0R*0Tl*0T2 

ITT*E«i 

GO  TO  30 
135  NS* 1 

rF(N.ca.ni  ns*i2 
r*i 

141  G(NSI*A(  11 

TF1G1NS1  .EQ.OI  GO  TO40 
N*NS 

IFtN.EQ.  11 1 GO  TO  SO 

IF ( N.EQ. 22) GO  TO  40 

NS*NS»1 

!*l»t 

GO  TO  140 

END 


I VJ*A|4I 


RI,R2,R3,R4  : Travel  time  in  seconds  from  zeroth  second  of  waves  reflected  one, two, three  or  four  times  from  the  bottom. 
01, G? UIO  : SetsMC  travel  ti»«es  ♦ruw  zerwth  second. 
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III 

COMPUTER  PROGRAMS  TO 
PLOT  DATA  (REFPLOTT , REDPLOT2) 
AND  PROGRAM  TO  COMPUTE 
STRUCTURE  (MULTLAY) 
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PP0G9AN  PfFPLOTT 

DIMENSION  XAXISCM  ,YAXlSI3>  , SHOT  ( tl  .AT  I Ul  .*(  <•! 

I NTfOEP  plOTFND,p$UP.PSDOHN,PFNUP, PENDONN 
PLOTCNO* -3  S PSUP*-1  1 “SOOWM*-? 

PFNUP*.1  ? PFNOO*N*2 
XAXlStl)  *<IMOTST AMC'E 
XAXISI2) *4H  2 SEC  P 
X A X IS  I 31 ■ SHE R INCH 
TAX  ISU»  t S 

Y AXISI2) *8HCC  ®E P I 
YAYISI  SI  *4HNCH 
9FAtM*,S)  XAXITI 4» , XAXlS(B) 

F FOPMATI  AS,AA«> 

PF  A0(9,9)T0TS 
o F OP  **AT <HX,FS.2) 

NX*TOTS/2.  *4  . 

IFf«X.L?.l J>  NX  a 1 3 

NXAXIS*40 

NY  A X IS*24 

CALL  PLOTINM-b.  ,6.  , 101 
CALL  PLOT  ( 24.,  0.  ,PLOf  ENOI 
CALL  POTATEXY 
30  10  1*1, NX 
x*r 

to  CALL  PLOTSYMBC  X,  J.  , .20,  13, 0.  ,PSCONM 
3 GO  dACK  TO  (4,-tl 

CALL  PL0TSYNPI4. ,-l.  . .2  4,  XAX  IS,0.,NXAX IS1 
CALL  PLOT!  0.  ,0.,p*NU»1 
00  ? 0 I*i,2F 
y«i 

20  CALL  PLOTSY*«l 0.,Y,.20, IT , 90. , FSCONN1 
C GO  BACK  TO  1-1.41 

CALL  PLOTSY»B<  -l  .,4.  , .24, Y AX  IS, 90.  .NTAXIS) 

JO  PFA0CS,%0»  ACE, ( AT1I1  , 1*1,111 
41  FO°**  AT  IA9,2F9.0,PF4.£l 
IFIEOF 141 IGO  T 0 BO 
IF<  ACf.EC.4H  I GO  *0  BO 

SHOT ( 1 ) * ACC 

ITINaATc  i) 

2«2f.-AT (2) • .00 1 •?. 

TS*AT|B| 

OIST=AT  f M >TS 
x *0 1 ST /?  • 

CALL  PLCTSY**fM  X,l.,.1  0,  S**OT  ,«<l  .,4| 

: T I **E  a T ( I ) * l SEC/INCH  ,OIST  *2SCC/INCM,AN0  x A N 0 v IN  TNChFS 
T T - 1 

oo  no  i*i,4 

PI  n -AT  I 1*71 
lFlPin.EQ.01  GO  TO  30 
100  1 1 * T 

SO  TO  30 
60  00  70  T*  1, 1 1 

IFIATII)  .E3.0.  I GO  T*  r 1 
y •*  A T f t ) *TS 

IFIY.GT.29.)  GO  TO  TO 
CALL  PlPTSYNBt  X, Y,.  t 0,1,0. ,PTur I 
r 3 CONTINUF 

CALL  PLOTSTNPfx,;,. « , |J,90.,PSUP» 

GO  TO  10 

43  XN«NX»4 

CALL  PLOT  f x N ,3  0 • ,Pl OT  rN0) 

FNP 
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PROGRAM  REOPLOT2 

DIMENSION  XAXlS<9).TAxISf5>  , SHOT  ( 1 » , A F ( j II  , 9 ( 41 
DIMENSION  OFPI 1 0 0)  «XD  f 1 00 1 

integer  plcten d,psup,ps do wn,penup*penoc*n 
PL  OTEND* -J  S PSUP*-1  t PSD0w*=-2 
«**rNUV»3  S PEN00WN*2  S DfPTw-0.  1 Nf=0 
KRITE  (61  ,31 

3 FORMAT  f t GIVE  REDUCING  VELOCITY,  F5.2*) 

READ  160* 4) REOVEL 

4 FORMAT (F5.  2) 

X AX  I S I ll*8WOISTANCE 
XAXISI2I*8M  1 SEC  P 
XAXISI3J *«M€R  I N CM 
TAXIS! 1> *8MTIME  .5 
TAX  IS f 2) *9HSEC  PER 
VAXISfM  *8H INCH 
VAXISC4I *OMT  - X / 

RE  AO  <5 ,0 ) XAXIS! 4)  ,XAXIS(5) 

8 FORMAT  f A A*  AR) 

PEAD!S,R)T0TS,V1 

9 FORMAT C15X,F5. 2, 5X,FS.JI 
NX*TOTS/ 1.T4. 

IF (NX.LT *13) MX *13 
MX  A X IS  *4  0 
NYA  X IS*4  0 

ENCODE (0,101, YAXTS (5) IREOVEL 
101  FORMAT (E5. 2*3X1 

CALL  PL0TINTI-6* ,6., 10) 

CALL  PLOT (29., 0* , PL C TEN 01 

call  rotatext 

00  10  Is 1 * NX 
X*I 

to  CALL  PLOTSTMR(X,0.*.20, 13*0. ,PSOOMNI 
C GO  RACK  TO  l%f-11 

CALL  PL0TST»**(4.*-1.  **25*XAXIS*0  .*NXAXIS> 

CALL  PLOT(0. ,0.*PENUP) 

00  20  1*1,27 
▼*I 

20  CALL  PLOTSTm®(0. ,T*. 20, 13*00. ,PSCONM) 

C GO  SACK  TO  l-Y**> 

CALL  PL0TSYM«(-l.tA.,.25,TAXIS  ,90., NY AX  IS) 

JO  p£AO(5,401ACE,(AT(X>  * I*  1,111 
40  FORMAT CAR.2FS.  0, or*. 21 
IF (EOF (5 ) I GO  T0  00 
IFCACE.EQ.5H  1 GO  TO  60 

SHOT ( 1 ) * ACE 
TTXPsATCl) 

NO*MOT1 

OEP < NOT*  27.-AT (2)*. 0 01*2. 

T S*  A T ( 6 1 
01ST*AT  < 7) ♦ TS 
RE0L’CC*0  1ST*  VI  7 RE  Of  EL 
X = OlST/t  . 

XO(ND) =X 

IFIX.GT.J7.)  GO  TO  55 
CALL  PLOTSY  m»| x , 1. *. 10, SHOT  *°0  .,91 
C TlMt  = fCI|«.S  SEC/IMCH.OTSTa?  SEC/ INCH  * ANO  X ANO  Y IN  INCHES 
55  T I*  1 

DO  100  1*1,4 
P f I • *A  T ( I ♦ 7 ) 

IFfRfD.EQ.OI  GO  TO  30 

too  n*i 

GO  TO  30 
60  00  70  !*l*lt 

IF  I AT  ( I ) • E Q • 0.  I GO  TO  30 
Y*AT ( I I »TS 
Y*(T-REOUCE) *?.A5. 

TF(V.GT.?7.)  GO  TO  70 

IF  (X.GT.37.)  GO  TO  70 

CALL  PLOTSYM8(X,T,.tO,3,0. ,PSUP) 

71  CONTINUE 
GO  >0  30 
00  oo  oo  1*1, NO 
X*XCfI) 

oepth=oep( n 

90  CALL  PLOTSYMPfx .DEPTH, .1,13, 90., PSUP1 
XNsNX»4 

CALL  PLOT(XN,JO. ,PLOTfNO» 

END 


C •' 

\ 

I ; 


I 
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PROGRAM  NUITIAT 

OlMENStON  Wf  ?<!!  , VI  20)  ,VAI  20!  ,VB4  20!  ,ALPH|?n)  ,EETA|  20!  , C4  201  , 

1 4 1201 ,01  201  ,T  »I(  200  ,T0T  120!  «HM  201  , MB  1 20!  ,0  A 1201  ,C0f  201  ,PC20I 

DIMENSION  TITLEIA) 

C IF  TNI  I NT  f 9CI.PT  TINES  APf  NOT  READ  IN,  I N$ERY  A BL  ANN 
C C A AO  WHERE  THtf  ARE  CALLEC  FOR  IN  iMt  CAT  A PECK.  THEY 
C WILL  THEN  PC  COMPUTED  BV  FORMULA  AFTER  ST  F°  421.  IF  THEY 

C ARE  READ  TN,  THEY  will  be  checked  FOR  ccns istfncy  within  10  PC 
402  WRITE  € f»  1*601! 


401  FORMAT  |#0«! 

C N* NUN BE  R OF  LAYFRS  OR  TRAVEL  T INF  SEGMENTS,  X*€NO-TO-ENC 
Z SPREAD  LENGTH 

®EAO  (5,405!  N,X,  (TITLE!  I)  ,IM  ,61 
IF<E0E(6))  CALL  EXIT 
400  FORMAT  < r*,F8. 0,6A8t 
IF  (N»  640,660,407 

402  REAP  (6,410!  l V A ( I!  , I - 1 , Nl 

610  FOR-AT  49FO.O) 

RE  AO  16,410!  4V8tlT,I  * 1,N! 

T A I 4 1 1 * .8 

°E  A P (6,4101  4 T A II  I)  , I*  2,  N) 

T 4 I 4 1 ! * .0 

read  <5,4io!  crercn  ,i*2,n> 

WRITE  (61,615!  Y,ITITLC4n  ,1*1,61 
415  FORMAT  (#  * , *°EVERSAL  DISTANCE  * # ,F 1 0 . 3 X , F A*, //! 

WRITE  461,420! 

420  FORNAT4 # INPUT  DATA#//#  LAYER  A Po  A RE  NT  APPARENT#, 

It  INTERCEPT  INTERCEPT  #/ 1 VELOCITIES  A vELO#, 

2FCITIES  B TINES  A TINES  B#l 

00  425  I*2,N 

f«0  * TA  1 4 1 ! ♦ XM1./V4II)  - 1./V0IIM 
IF  CTBI4TI!  622, 422,423 
422  TBI4IT  * TB8 

GO  TO  425 

421  TAENO  = TA  1 4 1!  ♦ XSVA(I) 

TREnO  * TBIIII  ♦ X/VP4II 
ERROR  * ABSE4TAENOFT9ENO  - 1.! 

IF  (ERROR  - .101  425,424,624 

424  WRITE  (61,14261  T 

1424  format  (#0APPARENT  VELOCITY  ANO  TINE  INTEPCEPT  OATA  ARE#, 

1#  INCONSISTENT  AT  LAYER  #/#  NUMBER  #,I2,#  ENO-TO-ENO  #• 

2 #TR  AVEL  TINES  DIFFER  °Y  MORE  THAN  10  PERCENT.#/! 

GO  TO  475 
475  CONTINUE 

WRITE  (6  1,  1426!  ( I , V A 4 1 ! , V°  (I!  , T A III!  , TB  I ( I » , I * 1 , * 1 

1425  F OR  NAT  4 # f , !6,r 1 3. 3,F  14.3,F12« 3, M 0.31 
VC1!  * 4VA41!  ♦ VB<t!l*.5 

00  57*  N * ?,N 
K * 1 

ALPHdt  * ASINF  4 V4  t ! /VB4N!  I 
BET  A 4 1 ! * ASINF4 V41I /VA(N1> 

IF  (*-?)  580,600,510 
*00  A 4 1 ! * 4ALPM41!  ♦ BET»(1M*.5 
W 4 2 ! * 4ALPM41I  - NETAIt!|*.5 
V42I  * V4  t! /SINE  CAT  1)  ! 

GO  TO  6*3 

510  A ( 1 ! * ALPM41!  - W ( 2! 

P 4 1 1 * BETA!  II  ♦ WC2I 
520  < * KM 

VV  * V(K)/V<K-lt 

PfK)  * ASINF  IV  V*STNE  ( A 4 K**  1 ! ! ! 


t 

ft 

m 


• J 


i 


1 


Q(Kt  * ASINF  fVtf*SINF<f  IK-1)  II 
IF  IK*1-"1  530 « 5*0 ,56  0 
550  UK)  * POO  - M IK*1I  ♦ MtK) 

R C K t * QIK)  ♦ WfK*l)  - NIK) 

ALP  H ( K)  a A ( Kl  ♦ tf(K»l) 

BET  AIK)  * 9IK)  - W(KM) 

GO  TO  5?  0 

550  Ann  = (piKi  ♦ Q«n*.5 
9(K)  X AIK) 

M « K » 1)  * M(K>  ♦ t P ( K ) - Q ( K ) ) ••  5 
ALPH(K)  X a no  ♦ M f K+  1 ) 

PETA(K)  * 9 00  - MIKM1 
V I K ♦ 1 ) * V(K)/SINF(A( *1  ) 

550  < K « K-t 

MHAxO. 

HH0  x 0. 

IFIK.E0. It  GO  TO  561 
00  560  T * 1,KK 

HM  X COSFCALPHf  in  ♦ COST  CPETAC  I)  I 

MH  x HM/ V( IT 

HHA  x HH A ♦ MM*MA|I> 

560  MM9  x HHB  ♦ MH*  H 0 C I ) 

561  9 * Vf  Kt /(COSFT ALPH(K) ) ♦ C OSF I9ET A (K t ) ) 

HACK)  = D*(TAIIK»ll  - HMAI 

m9CK)  * P*  f T*I ( K +1 ) - MM0) 

0 A 1 1 ) x HA(1) 

open  * nun 

0AIK)  x OA(K-t)  ♦ MA  ( Kt 
09151  * 08(5-1)  ♦ M8C51 
570  CONTINUE 

PT0Oxl90./3.  1615926 
00  5 AD  J x 2,M 
MO  M ( J ) x W(J)*PTOO 
49ITE  (61,620) 

620  F09  MAT  f * OCOMPOTE  0 ST® UCTUPE  tf/t  LATER  VELOCITY 

1 AT M IC5NF 9$  A THICKNESS  8 OIP  3EPTM  A CED  TM  **l 

r = i 

WPITE  (61,625)  I,  TII),  MA(I),  H0U),OA(I1,  Ct(H 
625  rOPHATU  t , 16,  F 1 2.  J , F It  . JfF  1 6.  3 , 9 X , F8 . 3,  F9  . J ♦ 

NN  x n-1 

MPT  TE  (61,636)  (I-  TO,  HUI),  Me(I),  M ( I ) ,C  A ( I)  , Cf  ( I)  • 
1 T = 2 * NN1 

633  FORMAT  ( t #,I6,F12,3,FU,3,F16.3,F9.3,F8.T,F9.3) 

MPITE  (6 1 * 635) N , VlNt,  W (N) 

635  c OP  M AT  it  f,I6,Fl2. 3,23X  ,F11.J) 

GO  TO  602 
660  CONTINUE 
END 


P 


i 

/■ 

i 


N : Number  of  layers  (equal  to  number  of  velocities). 

* : Ena  to  end  profile  length. 

riTlE(I)  : Line  identl flcatlon. 

1*1 N : Layer  mauler. 

VA(I)  ; Apparent  velocities  from  the  A-end  of  tne  profile. 

VB ( I ) : Apparent  velocities  from  the  8-end  of  the  profile. 

TA!(I)  : Intercept  times  for  layers  2 through  N at  A-end  of  the  profile. 

TB 1(1)  : Intercept  times  for  layers  2 through  N at  9-end  of  the  profile. 

All  distances  in  meters  or  kilometers. 

All  velocities  in  meters  per  seconds  or  kilometers  per  seconds. 

All  times  in  seconds. 

With  single-ended  line  : fB([)  • VA(I)  and  TB 1(1)  ■ TAJ([)  . 


OUTPUT  FORMAT  FOR  PROGRAM  MULT  LAY 
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IV 

INTERCEPT  TIMES  AND  APPARENT 
VELOCITIES  (INPUT  OF  MULTLAY) 
AND  COMPUTED  STRUCTURES 
(OUTPUT  OF  MULTLAY) 

FOR  ALL  LINES 


% 21 

.62  *GLt 

t.*% 

2.1% 

2.77 

*•0  2 

1 •%% 

2.1% 

2.  77 

9.02 

0.0%% 

0.  000 

0.680 

0.0%% 

0.060 

0.660 

REVERSAL  01  ST  AMCC  « 

21.620 

%GL1S 

INPUT 

LAVER 

DATA 

APPARENT 

APPARENT 

I NT  EXCEPT 

INTERCEPT 

VELOCITIES  A 

VELOCITIES 

9 TIMES  A 

TIMES  9 

1 

1 .%%# 

l.%%0 

0 

0 

2 

2.1*0 

2.  1*0 

• 0%% 

. 0%* 

3 

2.770 

2.  770 

. 000 

• 600 

% 

9.370 

9.  020 

.660 

• 0 60 

COMPUTED  STRUCTURE 

LAVER  VELOCITY  7NICKNCSS  A 

THICKNESS  9 

DIP  DEPTH 

1 

1.  *%• 

.0%I 

• 0*3 

• 1*3 

2 

2.  1%0 

. •%« 

• •%% 

0 .092 

j 

2.  770 

.%*• 

.9*6 

• 1.096 

% 

9.  026 

0 

DEPTH  P 
-f  4 3 
• 0 92 
1.090 


* 12. 77  «*(.* 

l.%%  2.1%  2.0% 

l«%%  ?.1%  2.0% 

0.031  0.230  I. 790 

0.031  0.230  0.7*0 


%.?€ 

9.2* 


REVERSAL  OISTAHCf  • 12.770  **01 2t 


rx.yr  otn 

LATEP 

APPARfN* 

APPARENT 

i*te *cr< 

VELOCITIES  A 

fELCCI^IES  f 

T 1 S 

1 

t.%%0 

l.**0 

0 

2 

2.1*0 

2.  1*0 

• 031 

3 

2.990 

2.  e9t 

• 239 

% 

9.269 

9.  2*9 

.7*0 

COPPUVCO 

i a ver 

ST»UCTu*e 

WIXTTY  T m ICKNC  S 5 t TMTC 

*.i€SS  9 

1 

l.**0 

. 9 JO 

• 039 

2 

2.  1*0 

.291 

• 296 

1 

2.  991 

• 6*0 

• •%• 

* 

9.2*0 

9 

intercept 

TT*ES  e 
0 

. 031 
.230 
.7*0 


OTP 


OCPTM  0 

• 030 

• 3?* 
1.17% 


OfPTH  » 
.030 
.3  2* 
1.17% 


5 

0 
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6 

36.61  SI 

9 TO  1 

•UZ*f 

1 .59 

3.35 

3.96 

1.01  J.777 

5.397 

1.59 

3.35 

3.9b 

S.  0 1 3*716 

5.359 

8.18 

9.  IS 

9 . 2% 

».k*  1.69 

0.  10 

0.19 

9.2  0 

9.63  1.29 

REVERSAL  01  ST  *NC€  * 

36.618 

91  6 TO  t M267S 

INPUT 

OAT  A 

layc® 

APPARENT 

APPARENT 

INTERCEPT 

intercept 

VELOCITIES  A 

VELOCITIES 

8 TINES  A 

TINES  8 

1 

1 .990 

1.998 

0 

0 

2 

2.2%  0 

3.398 

. too 

. 100 

3 

2. *60 

3.968 

.150 

. 150 

9 

3.310 

3.  818 

.300 

.23  0 

c 

3.737 

3.716 

• 95  8 

.930 

6 

9.387 

9.259 

1.900 

1.  290 

CONFUTED  STRUCTURE 

IA»ER 

VELOCITY  Tm 

ICKNESS  A 

thickness  r 

OTP  OEPYm  • 

OEPYm  8 

t 

l .690 

.899 

• 099 

.099 

• 099 

3.790 

.129 

. 125 

0 .317 

.217 

3 

3.980 

.037 

.027 

0 .299 

.299 

9 

3.  010 

• 57  8 

.537 

0 .833 

.771 

9 

3.  7? l 

2.18S 

1. 999 

• 118  3.008 

3.715 

6 

9.  333 

.629 

6 29.39  91 

9 *0  3 

*•395 

1.99  3.29 

3.56 

?.oo  3.898 

9.909 

1.99  3.39 

2.56 

2.99  1.867 

<••679 

3.10  0.16 

9.21 

8.56  1.375 

0.10  0.16 

0.20 

0.  58  1.  388 

RE Vc  RS  At  DISTANCE  * 

29.390 

•I  9 TO  3 -395S 

INPLT 

DATA 

LATE*? 

APPARENT 

APPARENT 

intercept 

INTERCEPT 

velocities  a 

VELOCITIES 

9 TINES  A 

TIPES  9 

1 

1 .990 

1.990 

9 

0 

3 

3.390 

2.290 

• 100 

. 100 

t 

2.560 

2.56  8 

. 260 

. 260 

9 

2.990 

3.  990 

• 200 

.200 

5 

3.899 

1.867 

.560 

.5  00 

6 

9.909 

9.679 

1.375 

1.080 

COMPUTED  structure 

later 

VELOCITY  TM 

ICKNESS  A 

YNICKNCSS  9 

DIP  DEPTH  A 

DEPTH 

1 

1.990 

* 093 

.093 

. 093 

.893 

3 

2.290 

.137 

.137 

0 .230 

.233 

3 

2.  560 

• 029 

.025 

0 .355 

.255 

9 

2.998 

.771 

• 630 

0 1.026 

• 885 

5 

3.  §92 

2.61  9 

1.6  83 

.276  3.995 

2.568 

6 

9.  789 

1 

.707 

* 

k 


» 


r 


ML. 
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5 10*5 

a.  arc 

WEST  76  9FAIF09T  NOA9A 

1640. 

1600. 

22J0. 

3929.  5690. 

1**0. 

1600. 

2230. 

3920.  5600. 

• 0 19 

.215 

.5  75 

t.  255 

.010 

.215 

.575 

1.255 

REVERSAL  CISTANCE  ■ 

10950. 

301  4tC 

WEST  75  PEAU*  OPT  P0A9AS 

INPUT 

DATA 

LATER 

APPARENT 

APPARENT 

INTERCEPT 

INTERCEPT 

VCL7C  IT  TES 

A VELOCITIES 

e tines  a 

TINES  t 

1 

1660.300 

1660.  001 

« 

0 

2 

1600. 109 

16  00.  000 

• 015 

.010 

3 

2230  .000 

2230. 000 

.215 

.215 

5 

3020. 000 

3020. 300 

.575 

.575 

5 

*640.900 

6600.  000 

1.255 

1.255 

CONPUTEO  STRUCTURE 

LATER 

VELOCITY  1 

mickncss  a 

THICKNESS  9 

DIP  oeptn  A 

CEPTH  P 

1 

1660. 009 

16.514 

16.510 

16.510 

16.510 

2 

1609.  990 

22 6. 406 

226. 006 

3 263.326 

26 J. 326 

3 

2230. 300 

515.010 

519.810 

0 763.162 

763.162 

5 

3920. 000 

1100.723 

1110. 723 

01063.065 

1063.065 

5 

6600. 009 

0 

r 2304*,. ARC  WEST  76  5EAUECRT  M099* 


1660. 

1600. 

2230. 

6260. 

5630.  *070. 

1660. 

1600. 

2230. 

29  30* 

6260. 

5*30.  7370. 

.010 

.215 

.670 

1.660 

1.925 

.310 

.215 

.670 

t.  MO 

1.660 

1.925 

Pf/*RSAt  DIST 

ance  ■ 

23066.0  00 

m 

WEST  76 

BEAUFORT  W0BR8S 

INPUT 

DATA 

LAYER 

APPARENT 

APPARENT 

INTERCEPT 

I NTERCEPT 

VELOCITIES 

• VELOCITIES 

B T I *CS  A 

TI«*ES  B 

1 

1660.000 

\64f. 900 

0 

0 

2 

1*00.300 

1609.  000 

.010 

• 810 

3 

2230.000 

2230. 008 

.215 

.215 

6 

2030.000 

2930. 000 

. 670 

. 670 

5 

6240.900 

6260. 000 

1.050 

t.050 

6 

5630.000 

5670. 009 

1.669 

1.660 

7 

7070. 330 

7 Of  0.  000 

1.925 

1.  925 

CONPUTED  structure 

LATER 

VELOCITY  T 

HICKNESS  A 

THICKNESS  9 

fll*  OERTm 

A CEPTH  ■» 

t 

1660. 000 

16.51  f 

16.518 

16.510 

16.510 

* 

1600. 000 

226. 006 

2 <6. 0 06 

0 263.326 

263.326 

3 

2230. 000 

391.099 

391.099 

0 635. 223 

635.223 

6 

2930. 000 

607.279 

067.279 

91522.501 

1522.501 

5 

6260.000 

906.128 

906.129 

02420. 6J9 

2628.629 

6 

5630. 000 

1750.796 

17*0. 706 

06179.615 

6179.615 

7 

7070.000 

0 

59 


7 2? 

.37  LI  *€ 

10-11 

12-13 

*ITn  6.0 

ANO  7.0 

L A YE°S  ' 

1.440 

1.650 

2.  *0  9 

3.0  44 

4.  056 

6.000 

6.406 

1.44Q 

t . 65  0 

2.509 

3.  090 

3.004 

5.000 

7.129 

0.009 

0.440 

0.675 

l.  130 

1.760 

1.805 

1.  008 

0.440 

0.  92  0 

1.360 

t.  760 

1.9  10 

REVERSAL  CISTANCE  * 28*070  LINE  10-11  12- 1J  KITH  t.O  AND  7.0  LAYERS  "1012S 


INPUT  04  T A 


LAYE  9 

APPARENT 

APPARENT 

INTERCEPT 

INTERCEPT 

VELOCITIES 

A VELOCITIES 

B TINES  A 

tines  p 

1 

1 .440 

1.440 

9 

0 

2 

1.65  0 

1.650 

. 00  9 

.009 

3 

2.509 

2.509 

• 440 

.440 

4 

3.044 

3.  090 

.675 

.020 

5 

4.856 

3.  884 

1.130 

1.  368 

6 

6.000 

6.  C 80 

1.  760 

t.  760 

7 

6.906 

7.129 

1. 80  5 

1.  830 

:o*puteo 

STRUCTURE 

LATER 

VELOCITY 

THICKNESS  a 

THICKNESS  B 

OIP  depth  a 

DEPTH  B 

1 

1.  440 

. 012 

.012 

.012 

.012 

2 

1.650 

.467 

.467 

0 .479 

.479 

3 

2.589 

.400 

. 716 

0 .879 

1.195 

% 

3.  067 

. 877 

. 966 

0.611  1.756 

2.161 

5 

3.  921 

1.085 

.362 

2.615  2.841 

2.523 

6 

5.  994 

-0.186 

.063 

Q.92T  2.855 

2.586 

7 

7.  306 

- 

2.  086 

6 78736. ARC  PEST  76  fEAUf"PT  «1617 

1460.  lSni.  1820.  2440.  3620.  6000. 

1440.  1600.  1 020.  24*0.  3620.  6000. 

.0  10  .100  .300  1.028  1.746 

.010  .101  .3*0  1.026  1.746 

PE  /5®SAL  CIST  A TCE  = 28736.000  ARC  KEST  76  PE  A cE  0 RT  -1617S 


INP'JT  34  T 4 

L A TE  8 APPARENT  APPARENT  I NTERCE*T  IiJTCPCEPT 

VELOCITIES  4 VELOCITIES  9 T’-fS  4 TINES  9 
1 1440.000  1440.  000  0 0 

? 1600.000  1600.  000  . 010  . 010 

3 1070.900  1020.  000  .100  .100 

4 2440.900  2440.000  .300  .300 

6 3620.000  3620.000  1.028  1.028 

6 6000.000  6080.  000  1.746  1.  746 

coh»uveo  STRUCTURE 

later  velocity  thickness  a thickness  b oip  ?cpth  a depth  r 

1 1440.  000  16.51  0 16.81  ! 16.518  16.810 

2 1608.  000  144.306  144.  306  9 160.  024  160.024 

3 1 920.  000  3 07.807  3 E7.807  0 460. E31  460.631 

4 2440.  000  907.846  507.  945  01376.  578  1376.575 

5 3620.  000  1226.193  1226.193  02602.76!  2692.760 

6 6000.  000  0 
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1 


6 26699 

.ARC  wf  *"T  7 S 

9f  A L F CRT 

1440. 

160  9. 

1 9 4 3. 

232t. 

3310. 

6660. 

! 440. 

1600. 

1 849. 

2323. 

3219. 

5860. 

.010 

.1  25 

.465 

1.  205 

2.311 

.313 

.125 

• 4 65 

1.205 

2.3  10 

reversal  distance  * 

26699.0  90 

A?C  REST 

76  9EAUF0RT  H1819S 

INPUT 

DATA 

LATER 

APP  ARENT 

APPARENT 

INTERCEPT 

INTERCEPT 

VELOCITIES 

A VELOCITIES 

8 TIMES  A 

TTPES  b 

1 

1440.000 

1440. 000 

0 

0 

2 

1600.000 

1600.  000 

.010 

. 010 

3 

1940.000 

1040. 000 

. 125 

. 125 

4 

2320.000 

2329.090 

.465 

.465 

5 

3310. 380 

33 1 1. 300 

1.205 

1.205 

6 

5660.000 

5660. 000 

2.  310 

2.310 

! 

COMPUTED  STRUCTURE 

LATER 

VELOC  ITT 

thickness  a 

thickness  9 

DTP  OF  PTH  A 

DEPTH  » 1 

i 

1440.000 

16.518 

16.518 

18.519 

16.518 

2 

1600. 090 

1 79.  367 

175.367 

0 196.  885 

195.885 

3 

1840.  3 0 C 

4?9. 958 

4 19. 958 

3 625.843 

625 .843 

4 

2320. 000 

975.136 

975. 136 

01600.979 

1631.979 

5 

3313. 090 

1762.815 

1 762.815 

33363.7R4 

3363.794 

6 

5660. 300 

0 

I 

6 22651 

. ARC 

NEST  76 

9EAUF0RT 

M?  32 1 

1 44  0. 

1600. 

1 890. 

2350. 

3 J90  • 

6680. 

1440. 

1600. 

1850. 

2350. 

3390. 

5688 

.010 

.125 

.490 

1.  105 

2.3  75 

• 013 

. 1 25 

.4  90 

i.  105 

2.175 

C€WF«>S*L  DIS TANCF  = 22651.000  ARC  mEST  76  ?*ALFORT  P2021S 


INPUT  .TATS 


LATER  apparent 

VELOCITIES  A 
1 1440 .30  9 

? 1600.000 

1 1890.310 

« 2360.300 

5 3390.090 

6 5690.000 
CONPUTEO  STRUCTURE 


APPARENT 
VEL JCITTES 

1 440. 000 
1600. 010 
l«90. 000 

2360. 000 
JT93. 000 

5690.000 


INTERCEPT 
9 TINES  A 
1 

. 019 
. 125 
.490 
1.106 
2.075 


intercept 
TIMES  ° 

9 

. 310 
. 126 
.490 
1.106 
2.  075 


LATER  VELOCITY 

1 1440.010 

2 1609.000 

3 1990.303 

4 2350.000 

5 3399.000 

6 6600.300 


THICKNESS  A 
16.51 8 
1 t5.535 
609.354 
740.227 
16  M.  993 


THTCK4ESS  8 

16.519 

1 15.535 
509. 354 
740.727 
1630. 983 


DIP  DEPTH 
16.518 
3 192.053 
3 691.487 
01432.  134 
0 3063.  1 17 


A DEPTH  9 
16.518 
182.053 
691.407 
1432.134 
3C63.117 


n 


. 
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6 44770 

. ARC 

wF  ST  76 

91  A l F 09 T 

•*22  024  0 

reverse 

1440. 

i6no. 

t 932. 

; 340. 

2932. 

4 79  1. 

1 *40. 

1600. 

1826. 

2316. 

2 780. 

6 034. 

• 0 10 

.1  76 

• 4 4 0 

. 010 

1.7  16 

.010 

. 1 00 

. 1 70 

.626 

2.  360 

RtVFRSAL  DISTANCE  • 41300.000  ARC  NEST  76  P^AIFCRT  P24A22RS  RtvrRSf 


INPUT  TATA 

LAYER  APPARENT 

VtlOCTTIFS  A 

1 1440.000 

2 tftSO.OQO 

3 19Y7.000 

4 JJ9«.000 

6 *.647.000 

CONPUTFO  STRUCTURE 

LAYER  ' VELOC TTY  THICKNESS  A THICKNESS  0 OIP  OtPTH  A DrPTM  4 

t 1660.000  16.  St*  16. MO  16.610  16.610 

? 1600.000  167. 8J9  282.007  o 166.167  299.604 

3 1 94  0 • 90  J 1 160.6  07  6 ° 7 . 1 76  -0.  1 72 1 324 . 764  1 *96.791 

6 3396.  660  1 069.669  .'032.60  0 . 03*3196.  4 13  3326.169 

5 6710.266  -0.166 


6 41000. ARC  WEST  76  Rf  • t'F  OPT  "24A22P  QTVERSf 

1 460.  IbOO.  1937.  JJ9S.  4*07. 

1460.  1600.  1966.  J 3 9 0 • 6 760. 

.010  .120  1.166  2.210 

.010  .216  1.179  2.290 


>?FVF9cAl  DISTANCE  * 46770.000  ART  wiST  76  RRAUFCRT  ■22A24PS  REVEST 


r * P t'T  TATA 

lAYF*  APPARENT 

VFIOCITIFS  A 

1 1440.000 

2 1600.090 

3 10  32.000 

4 234R.I00 

6 2932.000 

6 4 7,61.0  00 

CCPPUTFO  STRUCTURE 

LAYER  VELOCITY  THICKNESS  A THlC.NNFSS  R PIP  DEPTH  a OFPTH  « 

1 1640.000  16.6  1(1  16.610  16. M*  16.616 

2 1600.000  266.613  141.769  0 202.131  140.207 

3 1420.907  260.927  J4.027  .170  463.060  197.314 

4 2331.310  444.464  641.024  .4661096. M2  933.338 

6 2904.706  1100.006  2l«7.«0t  .3242190.690  3031.130 

6 4903.304  -1.036 


AP*A  R*"  NT 

I N T t RCF  ®T 

I NTF  RCFPT 

VFIOCITIFS  P 

T I RES  A 

TTRES  P 

1440. OHO 

0 

0 

1600.  300 

.010 

. 310 

102E. 000 

.176 

. 100 

2314.090 

.44  a 

. 170 

2 990,000 

. 01  0 

.4  24 

FO  3 4.  000 

t • 714 

2.  040 

APPARENT  INTERCEPT  I NT  f RCF  P T 


VELOCITIES  R 

1440. 000 
1600.  000 
1944.  001 
3390. 900 

6740. 000 


TIRES  A 
0 

• 010 
. 129 
1.166 
2.210 


T T RES  0 
0 

.010 
.214 
t.  1 74 
2.  290 


PROGRAM  TO  COMPUTE  THE  OFFSET 
DISTANCES  OF  THE  RAYS  (WAREFRA) 
AND  INPUT  AND  OUTPUT  DATA 


FOR  PROCESSED  LINES 
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PROGRAM  WAREFRA 

3 THIS  PROGRAM  COMPUTES  PLANE  DIPPING  LAYER  MODEL 
C RAfiANFTERS.SECTTONS  OF  L*TERS  ALONG  WHICH  RAY  IS  REFRACTED 
C ANO  THE  CRITICAL  DISTANCE.  THF  InRUT  DATA  REDUIRES 
C KNOWN  DR  APPARENT  LAYER  VELOCITIES  <VA»  ANO  EITHER 
C INTERCEPT  TIME  from  REFRACTED  ARRIVAL  OR  LAVER  THICKNESS 
C VERTICALLY  UNDER  THE  ORIGIN.  LAYER  DIP  IS  ASSUMED  1 FRO 
C IF  LEFT  RLANK.  IF  LAYER  THICKNESS  IS  GIVEN.  GIVEN 
C VELOCITY  IS  TAKEN  fC  The  TRUE  LATER  VELOCITY. 

DIMENSION  V f 10 ) • HOR  I10),W(10),VA(10).TIA(10),HA(10) 
DIMENSION  B « 1 0 > « ALPHA  (10  9*  RE  T A (10),  «p  <10>,  A (10) 
DIMENSION  XA(i8l  ,X8C  l 01  ,YA(tO)  ,Y8ll«l  ,KC(lO)  ,P(10> 
DIMENSION  5 ( 10 ) 

RT OO  ■ 1 AO • / 5. 14159265 
3T0R  * 1 . / RT  DO 
C REA 0 HOR  CARO 

to  READ  (5,  100)  N * vtll,  x,  HOR 
IFIE0FC5))  CALL  EXIT 

C N IS  NUMBER  *F  layers,  v « 19  IS  VELOCITY  OF  FIRST  LAYER, 

C X IS  LENGTH  OF  LINE,  HOR  IS  IDENTIFIER  MAX  LENGTH  OF  0 

130  FORMAT  ( IS • F5. 3 , F5 . 2 , 1 9 AS ) 

IF  (EOF  (1))  CALL  EXIT 
IF  (N  .IT.  10)  GO  TO  15 

WRITE  (61,  110)  N 

110  FORMAT  ( X-N  larger  THAN  DIMENSIONS  N x<  ,T5 ) 

CALL  EXIT 

1*5  K * 2 

VA  (IT  * V <1T 
W (1)  * 0.0 

3 REA 0 SUCCEEDING  CAROS 

3 IF  mA ( I ) IS  NOT  GIVEN,  COMPUTE  mqoEL  FrQn  t ia  INTERCEPT 
READ  15,  101)  VA  (?),  TTA  (?),  W (2),  HA  (1) 

101  FORMAT  (4F5.3) 

W (2)  * W (2)  • DT OR 

“ETA  (1)  * ASINF  (V  (l)  / VA  (2)) 

A (1)  * " (1)  * "ETA  (1)  ♦ W (2) 

ALPHA  (1)  * A (l)  ♦ w (2) 

V 12)  = V ( 1 ) / SIN  (A  (1  ) ) 

IFtNA(i)  .NC.0.0)  V f 2) *VA(  2) 

C FIND  out  if  thickness  given,  if  NOT,  CALCULATE  VELOCITY  ANO 

c angles. 

IF  (ha  (1)  .tc.  0.01  GO  TO  10 
IF  (TIA  (2)  .NE . 0.01  GO  TO  25 
22  WRITE  161,  111) 

111  format  (/-HA  ANO  TIA  .ED,  0.0/1 
CALL  EXIT 

Z COMPUTE  DEPTHS  TO  LASERS 

25  HA ( l)*(T  TA(21*Vfl) )/ l COS ( AL “HA ( 1 ) )*COS (PET  A(1  ) | ) 

30  HR  (t)  1 H*  111  . x • f ANF I W (2)) 

C COMPUTE  CRITICAL  DISTANCES  ANO  ENO  POINTS  IN  LATER 
c where  ray  travels  AS  REFRACTED  WAVE 

XA  (1)  x ha  (1)  • TANF(ALPMA  (l)) 

1 / (t.l  ♦ TANF(W  (2)1  • TAAFIALPHA  (til) 

YA  (i)  x hA  (1)  - XA  (11  • TANFtW  (2)1 

XC(1)*CHA(1)-XA(1)»VANF(V(2)))*TANF(B(1)MVA(11 
X"  (1)  * MB  (1)  • TANFIPETA  (lM 
t / (1.3  - TANF(BETA  (in  • TANF  ( M (2)1) 

TB  (1)  X HR  (1)  » / B (t)  * T AN  F ( W (21) 

XN  (1)  * X - XB  ID 
40  K * K ♦ I 
3 DO  LAYERS  3 THROUGH  N 

pc  A C 15,  10M  VA  (K),  TTA  (K),  W (<l,  HA  (K-l) 

3 COMPUTE  ANGLES 

W (K)  * W (K)  • OTOR 

"ETA  (1)  x ASINF  (V  (1)  / /A  (Kll 

“(1»  * BETA  (1)  * w (21 

I = l 

KM|  x < - i 

50  I * l ♦ l 

0 (I)  * ASINF  (V  (IT  / V (1-1)  • SIM  IB  II-t))T 
BETA  (!)  x 0 (I)  - w (I) 

B (I)  * BETA  (I)  #•  w IT  M ) 

IF  (I  .LT.  KMt ) GO  TO  50 
A CK-I)  * B IK  - 1 ) 

ALPHA  (K-1T  * A (K-l)  ♦ W (K) 

V I K)  * V (K-l ) / SIN  (A  (K-l) ) 

IF  ( HACK- 11  .NC.O,  01  vikt^yaikt 

P (K-l)  x A (K-t)  - W (K-t)  ♦ W CK) 

1 * K - 1 

70  T x T - t 

A (I)  * ASINF  (v  (I)  / V (IM)  • SIN  (P  (IMM) 

AL°HA  (IT  * A (I)  » M 1 I PI 9 
P (I)  x a (I)  - w (D  • W ( t ♦ 1 9 
IF  (I  .NE.  1)  GO  TO  ?0 
KM2  x < - 2 

TEMPI  x TEMP?  » TEMPj  * 0.0 
00  00  I X 1,  K M? 

TEMPI  * TEM“t  ♦ HA  (I)  / f <11  • (COS  (AL®HA  (ID 
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t ♦ COS  (BETA  I II  >) 

IMP*  a TEMP?  ♦ Mi  (II 
TEMPI  * TE  MP3  * H0  III 
80  CONTINUE 

IF  (Hi  (K-l)  .NE.  0.0)  GO  TO  85 
C FJKO  OUT  IF  THICKNESS  GIVEN.  IF  NCT  COMPUTE  VcLOC  ITT  iNC 
C MORE  INGLES 

IF  (TIi  (<)  .EO.  3.  0)  GO  TO  2* 

C COMPUTE  DEPTH  TO  LAYER 

HI  (K-t)  a (TIi  m - TEMPI!  • V (K-l) 

1 / (COS  (ALPHA  (K-l))  ♦ COS  (PET I <K-1))I 

C COMPUTE  01  ST INC  FS  IN  LATER  WHERE  WAVE  IS  REFRACTEO 
85  TEMP2  * TEMP?  ♦ HI  (K-l) 

HP  (K-l)  * TFMP?  - TEPP3  - X • TANFfW  (K)| 

TEMP  3 a TEMP  3 ♦ H0  (K-l) 

TF^PXI  s HA  (l)  • TANFfALPHI  (1)1 

l / (1.0  ♦ TA4F(W  (2))  • TINFIALPHI  (1)1) 

TE-PT0  a H0  (1)  • TANFTBETI  (11) 

1 / (1.0  - T INF ( N (2)1  • T A NF ( BET  A (1)11 

00  1 15  1 = 2*  KH1 

TEMPXI  « TEHPKA  ♦ (HA  (X)  • TINFIALPHI  ( I)  1 ♦ TE*PXA 

1 • T ANF(  AL  PH  A (I))  • ( T INF  ( W (1)1  - T INF  ( H (I*l))>| 

2 / (1.0  ♦ T INF  ( ALPHA  (1)1  • TINFIW  (IM))) 

TEMPX8  * TEMPXB  ♦ (HP  (I)  • TANF (BETA  (I))  - TEmPX? 

t • TINF(PETA  (1)1  • ( T ANF ( W CM  - T A NF  ( W (1*11)11 

2 / (1.0  - TINF(PETA  (I))  • TINFIW  (Ift)ll 

115  CONTINUE 

XI  (K-l)  a TEMPXI 
X 9 (K-l)  a X - TFMPXB 

T I TK-1)  a TEMP?  - XA  (K-l)  • TANF(W  ( K)  ) 

TP  (K-l)  a TEMP3  ♦ TEPPXP  • TAKFIW  (Kl ) 

1 a K 

DELTAXC  a TE^XA 
123  I = I - 1 

OELTAXC  a OELTAXC  ♦ (HA  (X)  • T ANF ( PET  A (1)1  ♦ OELMXC 
t • TANFtBETI  (II)  • (TAAE(M(I)I  - TANF  ( W (IM)>  * 

2 / (1.0  - TANF  ( w (I))  • TINFOETI  (If)) 

IF  (I  .GT.  2)  GO  TO  120 

XC  (K-t)  a OELTAXC  ♦ (HI  (II  - OELTAXC  • T A NF ( W (2»»l 
1 • TANF(BETA  (in 

C GO  BACK  FOP  ANOTHER  LATER 
IF  (K  .LT.  N)  GC  TO  40 
C CONVERT  ANGLES  TO  DEGREES 
OO  130  I a I,* 

«(!)*«  (I)  » PTOO 
130  CONTINUE 
C PRINT  OUT  RESULTS 

WRITE  (61.  102)  MOR,  4,  X 

102  FORMAT  ( 101. 10  Aft./  . * M * * • 12. A SPPEAC  a *, F6.lt 
WRITE  (61.  103) 

103  FORMAT  I *0  N APPARENT  DIP  LAYER  THICKNESS#./. 

1 < TELOCITY  AT  ORIGIN #1 

N Ml  - N - 1 

00  140  I a 1,  NM1 

WRITE  (61.  104)  I,  VA  (II,  W (I),  hi  (I* 

134  fo**at  ( I3.FR.  J.FP.J.6X.F6.3) 

140  CONTINUE 

WRITE  (61,  1P4)  N.  VA  (N).  W (Nl 
WRITE  (61,  135) 

105  FORMAT  (#0  M DEPTH  A THICK  A VELOCITY  THICK  0 # 

1 #OEPTH  B* I 
SMA  a 0.0 
SHM  * 0.0 
00  150  I a 1,  MM1 
SHA  a $H  A ♦ HI  (I) 

SHP  a SMB  ♦ HR  (!) 

WRITE  (61,  106)  I,  SHI.  HI  (I),  V 1 1>  « HR  (I).  SHB 
lie  FORMAT  ( 13. F8.3,FR.3,F10.3,F10.J.F<*.3) 

150  CONTINUE 

WRITE  (61.  10F)  N,  V (N) 

187  FORMAT  ( 1 3 , 1 8X  , F 8.  2) 

WRITE  (61,  108)  X 

188  FORMAT  (# OLIYER  RAY  TO  LITER  CRITICAL  CIST  RAY  FROM  #, 

1 #L A TER#, /, t X Y X#, 

2*  T#./,<  8.00  C. 08  #, 

3 FT . 2 , # 8. 00  # I 

OO  168  X a l,  NMl 
IPl  a I ♦ i 

WRITE  (61.  109)  T,  XA  (I),  TA  (II,  XC  (II,  XR  (II,  YB  f I I 
10«  FORMAT  (T3«FP.3,F8.3,Ftl.3,Fii#3,F9.38 
160  CONTINUE 
GO  TO  10 
END 
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INPUT  FORMAT  FOR  PR06RAN  MARE  FRA 


54 

\n 

X 

TlE(I) 

wuij 

TAI(1| 

— 

1 * 4 1 

M(2J 

f A I { 2 ) 

M(2) 

■on 

■BH 

• • • ' 

1 i 

M(N) 

IKJlfll 

■ • ‘ ‘ ■ - . . v ^ 

N 

: Number  of  layers  (equal  to  number  of 

velocities) 

VI 

Mater  velocity. 

X 

End  te  end  st  'ad  length. 

title; i) 

: line  identification. 

I-l N 

: layer  number. 

VA(  I ) 

: Apparent  velocities  at  A-end  or  true 

velocities. 

TAIU) 

: Intercept  times  at  A-end. 

M(U 

-ayer  dips  in  degrees. 

hA(  I ) 

layer  thicknesses  at  A-end. 

All  times  In  seconds. 

All  distances  In  kilometers. 

All  velocities  In  kilometers  per  seconds. 

'he  orogram  needs  Intercept  times  or  layer  thicknesses. 

The  layer  dip  Is  assumed  zero  if  left  blank. 

If  the  layer  thickness  is  given  , the  given  velocity  is  taken  to  be  the  true  layer  velocity. 


-ITLEil) 

» 

I«l 

VA(  I ) 

4(1) 

HA(I) 

HB(  I) 

V(I) 

SHA(I) 

SH8 ( I ) 
<A(I) 

XB<  I ) 

*CU) 

«U) 

»B(I) 


line  identification. 

NumDer  of  layers  Including  the  water  layer. 

End  to  end  soreao  length. 

Layer  number. 

Apparent  velocities  at  the  A-end  of  the  profile  or  troe  velocities, 
layer  dips  in  degrees. 

Layer  tnicknesses  at  A-end. 

Layer  thicknesses  at  S-end. 

True  velocities. 

Layer  depths  at  A-end. 

Layer  depths  at  B-end. 

Offset  distances  from  A-end. 

Offset  distances  from  S-end. 

Critical  distances. 

Layer  depths  at  the  offset  distances  at.  A-end. 

Layer  depths  at  the  offset  distances  at  8-end. 


All  distances  In  kilometers. 

All  velocities  In  kilometers  per  seconds. 


61.66  31.6?  NGL1 
3.136  .044 
3.76?  .090 
5 . 0 1 T .6?# 

WGL1S 

N * 4 SPREAD  • 31.6 

N APPARENT  DIP  LATER  THICKNESS 

VELOCITY  AT  ORIGIN 

1 1.660  0 ■ .06J 

3 3.  136  0 .069 

3 3.  769  0 .°66 

6 S. 017  a 

N OEPTH  4 THICK  A TELOCITY  THICK  0 DEPTH  9 
1 .06  3 • 065  t.660  .063  .063 

> .093  .069  3.136  .069  .0°3 

3 1.097  .966  3.76?  .966  1.99? 

4 9.03 

ESTE»  PAT  TO  LAYEP  CRITICAL  OIST  PAT  PROP  LATE* 
XT  XT 

0.  00  0.  00  31.63  0.00 

1 .039  .063  .070  31.9*1  .043 

3 . 009  .093  .171  31.939  .0«3 

3 .679  1.097  1.390  30.969  1.097 


61.66  13.77  MGL2 
2.136  .031 
3.996  .330 
9.399  .760 

NGL2S 

NT  6 SP»fAO  * 13.9 

N APPARENT  01*  LATER  THICKNESS 

velocity  at  origin 

1 1.660  ? .030 

3 3.  136  0 . 396 

3 2.996  0 . *93 

6 9.399  9 

N DEPTH  4 THICK  4 VELOCITY  THICK  B DEPTH  9 

1 .030  .030  1.660  .030  .038 

2 .339  .396  3.136  ,396  .329 

3 1.172  .*93  2.996  .893  1.17T 

4 9.26 

later  rat  to  later  critical  oist  rat  prop  iatfp 

XT  XT 

0.  00  0.00  13.77  0.  00 

t .030  .030  .099  13.763  .030 

3 . 316  . 339  . 633  13.696  . 339 

3 .790  1.177  1 .661  13.060  1.177 
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61.44  26.61  3T  b TO  1 W1267 

2.249  . 130 
2.461  .150 
3.013  .200 

3.721  .Ilf*  .578 

5.322  . 6252. 183 


6 T3  t N1267S 

N * 

6 SPOEA0  - 

26.6 

N 

apparent 

DIP 

LATER  THICKNESS 

VELOCITY 

AT  ORIGIN 

1 

1.440 

3 

. 094 

2 

2.249 

0 

.122 

3 

2.463 

9 

. 028 

i* 

3.  013 

0 

.578 

5 

3.721 

.116 

2.  183 

6 

5.322 

• 62% 

N 

PE**TH  A 

THICK  a 

velocity  thick  e 

OEPTH  9 

1 

• 094 

• 094 

1. 440 

.094 

094 

2 

• 216 

. ir> 

2.249 

.122 

216 

3 

.244 

. 029 

2.463 

.02! 

2\4 

4 

.822 

.578 

3.013 

.524 

769 

S 

3.  005 

2.  183 

3.721 

1 .947 

2 

716 

6 

5.32 

LAYER  ° A Y TO 

LA  YEP 

critical  dist 

RAY 

LAYER 

X 

Y 

X 

Y 

0.  00 

0.00 

26.61 

0.  00 

1 

.078 

.094 

.156 

26 

.532 

.094 

2 

• 342 

• 216 

• 6 83 

26 

.266 

• 216 

3 

. 228 

.244 

.457 

26 

.38  2 

.244 

4 

.960 

. 920 

1 .912 

25 

.727 

.770 

5 

2.  716 

2.976 

5.323 

24 

• 219 

2.  742 

61 

.44  29.39  61  4 TO  3 

*345 

2.209 

• 100 

2.556 

• 168 

2.993 

. 210 

3.882 

.276  .771 

4. 794 

1.7072.419 

91 

4 TO  3 K345S 

M * 

6 SPREAD  — 29.4 

N 

APPARENT 

otp  l *te r thickness 

VELOCITY 

AT  ORIGIN 

1 

1 • 440 

0 .093 

y 

2.289 

0 . 138 

3 

2.556 

3 .023 

4 

2-  993 

0 . 77 1 

5 

3.882 

.2/4  2.419 

6 

4.  794 

1.707 

M 

3E°TH  A 

THICK  A VELOCITY  THICK  8 

OE»TH  9 

1 

. 393 

.003  1.66  9 

• 093 

.09  3 

2 

.231 

• ua  2.28° 

.138 

.231 

3 

.254 

.323  2.336 

• 0 <3 

.254 

4 

1.  025 

.271  2.833 

.629 

. 863 

5 

3.444 

2.613  3 • 882 

1.6  85 

2.568 

6 

6.78 

LATER  pay  TO 

LJYEV  CSITIC8L  OrST 

RAT 

PROP  LATER 

» 

T 

X Y 

0.00 

9. 09 

23.  33  0.09 

l 

.075 

.133  .139 

29 

.315  .393 

2 

.341 

.231  .682 

29 

.049  .231 

3 

.263 

.236  .366 

29 

• 13  7 .254 

4 

1.  134 

1.920  2.131 

26 

•464  .888 

5 

4.4J6 

3.312  8.336 

26 

•291  2.661 

70 


51.99  10 

.Ob  M8A9A 

1 .blf  .010 

.015 

2.230  .21* 

.227 

3.020  .576 

.520 

<*.bfl  (1.259 

1.  101 

**A9AS 

N a 5 

sp®c&o  * 10.9 

N 

AP® A RENT 

9 IP 

LAYER  THICKNESS 

VELOCIT  y 

or  origin 

t 

1.  ‘•^O 

0 

• 016 

2 

1.608 

0 

.222 

3 

2.  230 

3 

.520 

i» 

3.  020 

0 

1.  101 

5 

5.600 

9 

N 

DEPTH  A 

3*1  CK  A 

VELOCITY  THICK  P 

DEPTH  9 

1 

• 016 

.016 

1.590 

.0  16 

• 01  6 

2 

• 293 

.222 

1.600 

.222 

.293 

3 

• 263 

• 520 

2.230 

.520 

.763 

% 

1. 969 

1.  101 

3.020 

1.101 

1.669 

5 

9.69 

u«»f»  9*r  ro 

LAYER 

CRITICAL  OIST 

rat 

FR  0 * L * YfR 

X 

T 

X Y 

0.00 

0.  09 

19.  96  0.69 

1 

.033 

.016 

• 366 

10 

.527  .016 

2 

.292 

• 291 

.9  95 

10 

.613  .293 

3 

. 220 

. 263 

1 .990 

10 

.190  .?63 

% 

1.300 

1. 969 

2.599 

9 

• 560  1.969 

M.?*  ?3. 

95  W9890 

l.*OC  .(to 

• 016 

2. 23i  .Ft? 

. 222 

i.sjt  .?r a 

.392 

?.?Wl.9?0 

• 907 

?.6JCl.<.<»Q 

.906 

7.0MI  .121 

l.  751 

kobrrs 

N * ? SPREAQ  a ?3.9 


N 

A P 9 4 RENT 

DIP 

14YF»  r«IC«'*FSS 

VELOCITY 

AT  ORIGIN 

1 

1.  590 

0 

. 016 

2 

1.60Q 

9 

.227 

3 

2.230 

9 

.39? 

9 

2.930 

9 

• 087 

5 

9.  290 

0 

. 906 

6 

5.630 

1 

1. 751 

2 

7.  070 

3 

N 

3f ®T  H A 

Thick  a 

VELOCITY  Thick  f 

DEPTH  B 

1 

• 016 

• 016 

1.  590 

. 0 16 

.016 

? 

• 293 

• 227 

1.600 

.22  7 

.293 

3 

.636 

. 3 92 

2.230 

.397 

.635 

9 

1.522 

. 987 

2.  930 

• 897 

1.52? 

5 

2.920 

.906 

9.290 

.906 

2.928 

b 

9.179 

1.751 

5.630 

1 .751 

9.179 

7 

7.07 

LAYER  RAY  t0 

LAYER 

CRITICAL  OIST 

RAY  LAYER 

X 

V 

3 Y 

9.  00 

0. 01 

21 

• 95  0.  00 

1 

.033 

.016 

• 0 66 

23. 

917  .016 

2 

.297 

.293 

.9«6 

23. 

6t I .293 

3 

. 667 

.635 

1.333 

23. 

183  .635 

9 

1.136 

1.6  22 

2.271 

22. 

719  1.522 

9 

1.  793 

2.  920 

3 .587 

22. 

05?  2 .9  ? X 

6 

3.550 

9.179 

7.116 

2t. 

29?  9.170 

71 


7 1.44  33.6  *41011 


1.65#  .030 

2.  J2C  . 365 

3. Q1f  .670 
4. 1 • C 1 . 140 
5.6101.545 
6.9t3l.805 


N 

APPARENT 

DIP 

LAYER  THICKNESS 

VELOCITY 

ay  origin 

t 

1.440 

n 

• 044 

2 

1.650 

8 

.372 

3 

2.320 

0 

.42  9 

4 

3.  0J0 

0 

. 821 

4 

4.  100 

0 

. 787 

6 

4.610 

1 

.762 

7 

6.  910 

9 

N 

Of  **T  H A 

THICK  A 

VELOCITY  THICK  8 

DEPTH  9 

1 

• 044 

• 044 

1. 44  0 

0 44 

.044 

2 

• ••  16 

.37? 

1 .65  0 

372 

.416 

3 

• 845 

.429 

2.320 

429 

.845 

4 

1.666 

.421 

3.030 

921 

1.666 

5 

2.453 

. 707 

4. 100 

797 

2.453 

6 

3.2  35 

.782 

5.610 

7 82 

3.235 

7 

6.41 

LAYER  PAY  TO 

t a yer 

CRITICAL  OIST 

RAY 

FRO*  LAYER 

X 

Y 

X Y 

0.  00 

3.00 

3 0.  6 0 0.  00 

1 

. C79 

.044 

.157 

30 

.521  .044 

2 

.411 

• 4 16 

.822 

30 

•189  .416 

3 

.776 

• 845 

1.55? 

29 

• 824  • 8 AS 

4 

1.375 

1.666 

2.750 

29 

•225  1.666 

5 

1.691 

2.453 

*■.381 

28 

.909  2.453 

6 

2.322 

3.235 

4.644 

28 

.279  3.235 

*12tJS 
H * 7 


7 SP?E  AO 
»P£>I  RENT 
VELOCITY 
1.  442 
1.650 
2.540 
1.460 
6.  000 
5.350 
7.140 


71.44  32.4  R1213 

1.65C  .030 

.044 

2.542  .450 

.433 

3.460  .960 

.827 

4.0001.240 

. 965 

5.9801. 755 

.519 

7.1401.940 

.94? 

1 * 3?. 4 

DIP  layer 

THICKNESS 

Af 

ORIGIN 

DEPTH  a 

THICK  A 

velocity 

Thick  9 

DEPTH  4 

• 344 

.0  44 

1. 44  0 

.04-4 

.044 

.477 

.433 

1.65  0 

.433 

. 4 77 

1.304 

.827 

2.548 

• 8 27 

1.304 

2.164 

.8  65 

3.462 

.865 

?.16« 

2.688 

.514 

4.800 

.519 

?•  688 

2.730 

.042 

5.990 

.04? 

2.  738 

R RAY  TO 

layer 

7.  14 

CRITICAL 

OIST  RAY 

FRO*  LAYER 

X 

0.  00 

f 

9.00 

X Y 

32.40  0.  00 

.079 

. 044 

.157 

J2 

•321  .044 

• 400 

.477 

.801 

J< 

• 000  .4  77 

1.144 

1.304 

2.2  98 

11 

.291  1.304 

2.  384 

2.  164 

4.768 

10 

.016  2.164 

1.604 

? • 688 

3.208 

?0 

.796  2.688 

1.  321 

2.  730 

2.643 

11 

.0/4  2.730 

72 


61. W4  29.74  M lb  1 7 


1 .60  0 .31  0 • 01S 

1.420  .100  .144 

?.640  . 1*0  . 304 

I.62C1.025  .904 

6.0811. 7 49  i. 226 


Hlbir s 

N * 6 SP®EAO  a 29.  7 


N 

APPA  RENT 

3IP 

LAYER  THICKNESS 

VELOCITY 

AT  ORIGIN 

1 

1.440 

0 

• 016 

2 

1.600 

0 

.164 

3 

1.029 

0 

• 309 

4 

2.440 

0 

. 908 

5 

3.620 

0 

1.226 

6 

6.  080 

0 

N 

OFPTH  A 

THTCK  a 

VELOCITY  THICK  6 DEPTH  8 

1 

• 016 

.016 

1 • 44  0 

016  .016 

2 

. 160 

.166 

1.600 

144  .160 

5 

• 468 

• 3 08 

1.020 

3 C 0 .466 

4 

1.376 

.908 

2.46  0 

900  1.376 

S 

2.692 

1.225 

3.620  1. 

226  2.602 

4 

6 . 08 

LATER  9 AT  TO 

LAYER 

CRITICAL  OIST 

RAY  PRO  7 lAtF» 

X 

Y 

X Y 

0.  00 

0.00 

29.74  0.  00 

1 

. 033 

.316 

.066 

29.707  .016 

2 

.286 

. 169 

.5  73 

29.656  .160 

5 

.402 

.668 

.963 

29.259  .669 

4 

1.986 

t.  376 

2.171 

20.654  1.376 

5 

1.  446 

2.602 

2.0*2 

20.294  2.602 

61.44  2%. 70  M16t  9 
1 .409  .910  .916 

1.848  .129  .170 

2.328  .465  .430 

3.3101.205  .975 

5.8602. Jl I 1. 753 


-1019$ 


N * 

6 SPOCAQ  * 

26.7 

N 

APPARENT 

OIP 

layer  THICKNESS 

VELOCITY 

at  ORIGIN 

1 

1.460 

9 

. 016 

2 

1.  600 

9 

. 179 

\ 

1 . 840 

0 

.630 

4 

2.329 

3 

.976 

6 

3.310 

9 

1.763 

6 

5.  660 

9 

N 

DEPTH  A 

THICK  A 

VELOCITT  IMICt  9 

OE»TM  4 

1 

.916 

• 716 

1.46  0 

016 

.01  6 

2 

. 1 °5 

• 179 

1.600 

179 

.195 

3 

• 625 

• 630 

1.060 

439 

.625 

6 

1.6  no 

.975 

2.320 

975 

1.600 

6 

3.  363 

1.  763 

3.310  1. 

763 

3i363 

6 

5.65 

lay 

ER  RAY  TO 

layer 

CRIT  ICAL  OTSY 

RAT 

PRO*  late® 

t 

Y 

X Y 

3.00 

0.00 

26.70  9.  03 

t 

.033 

.016 

.0  66 

26 

• 667  .0  16 

7 

. 335 

. l<* 

.671 

26 

.365  • l ®5 

x 

. 743 

.625 

1 .486 

25 

.997  .625 

6 

1.352 

t .609 

2.70% 

25 

•340  1.600 

5 

t.  °14 

5.363 

3.020 

26 

.706  5.363 

61.44  22. 65  *2921 


1.60C.QIQ  .016 
1,590  ,125  .165 
2.35C  .590  .*39 
3.3901. 105  .751 
5.6802.075  1.631 


M2021S 

N 3 6 SPREAD  * 22.6 


N 

APPA  RENT 

DIP 

LATER  THICKNESS 

VELOCITY 

at  origin 

1 

1.560 

0 

• 016 

2 

1.600 

0 

.165 

3 

1.690 

0 

.509 

6 

2.350 

0 

.761 

5 

3.  Ja0 

5 

1.  631 

6 

5.660 

0 

N 

de»tm  A 

THICK  A 

VELOCITY  THICK  6 

0€*TH  9 

1 

• 016 

.016 

1.460 

.0  16 

• 01  6 

2 

. 181 

.165 

1.600 

.165 

.181 

3 

.690 

.5  09 

1.  890 

.5  C9 

.690 

w 

1.431 

.761 

2. 350 

.741 

1.431 

5 

3.  062 

1.631 

3.  J90 

1.631 

'.062 

6 

5.66 

LAYER  RAY  T r 

LAYP» 

CRITICAL  01 S T 

RAY 

PR  0 * LAYfR 

X 

Y 

X Y 

3.00 

3.00 

22.65  0.00 

1 

. 033 

.016 

• 0 *6 

22 

.617  .016 

2 

. 2 81 

. 1*1 

.562 

22 

.369  .171 

3 

. 655 

.6  90 

1.709 

21 

.795  .650 

6 

1.150 

1.631 

2.301 

21 

•500  1.631 

S 

1.782 

3. 062 

3.566 

20 

• 666  S. 0 6? 

€1 

.66  65.77 

N2?A?%6 

1.60C 

• 010 

• 016 

1.829 

♦ 0.17 

• 266 

2.331 

♦ 0.65 

• 201 

2.916 

♦ 0.32 

• 535 

5.903 

-1.031 

• 100 

M22A260S 

H * 6 SPREAD  * 4*.f 


N 

APPARENT 

OIP 

LAYER  THICKNESS 

VELOCITY 

AT  ORIGIN 

1 

1.660 

9 

.016 

2 

1.  600 

6 

• 266 

3 

1.  829 

. 170 

.281 

6 

2.33! 

.450 

.535 

5 

2.906 

. S20 

1.160 

6 

5.903 

-1.0  30 

N 

DEPTH  A 

THICK  A 

VELOCITY  Thick  9 

DEPTH  8 

1 

• 016 

• 016 

1.460 

.0  16 

• 016 

2 

#2  82 

.266 

i.eoo 

• 130 

• 166 

3 

.563 

• 2 81 

1.629 

.157 

.204 

6 

1.090 

.535 

2.331 

.639 

• 64? 

5 

2.198 

1.  160 

2.90  6 2 

.179 

3.021 

6 

5.50 

l*»e*  *»r  ro 

LAYFR 

CRITICAL  DIST 

RAY 

FROI*  LAYER 

t 

Y 

x y 

9.  09 

0.00 

65.77  0.  00 

l 

.033 

• 616 

.0  66 

45 

.737  .016 

2 

• SOS 

• 26t 

1.003 

65 

.513  .147 

3 

.633 

# 558 

1.267 

65 

.500  .205 

6 

1.  162 

1.092 

2.303 

66 

.752  .868 

5 

.967 

2.215 

2.014 

66 

.201  2.993 

II 


74 


51.46  HI. 30  424A??q 

1.600  .310  .01S 

1.941  -0.17  ,146 

3.316  *0.031.160 

6.716  -0.181.879 


R24A22*  s 

N « 6 SP®EA0  * i»i.o 


N 

APPARENT 

OIP 

LAYER  THICKNESS 

VELOCITY 

AT  ORIGIN 

1 

1.440 

0 

. 016 

2 

1.600 

0 

• 148 

3 

t.961 

-0.170 

1.  160 

4 

3.  396 

. "30 

1.  870 

6 

6.718 

-0.1*0 

N 

DEPTH  A 

THICK  A 

VELOCITY  THICK  8 

DEPTH  3 

1 

• 0 16 

.016 

1 . 4%  3 

0 16 

.316 

2 

. 1 64 

.148 

1.600 

273 

.286 

3 

l .324 

1.160 

1.941  1. 

0 17 

1.3P3 

4 

3.  1 ° 4 

1.  870 

3.396  2. 

320 

3.  323 

«; 

6.7? 

LAYER  PAT  TO 

LAYER 

CRITICAL  31 S T 

PA  Y 

PROP  l»YE® 

X 

Y 

1 Y 

3.  00 

3.  03 

4 1.30  0.30 

1 

. 03J 

.016 

.066 

40 

.<*67  .316 

2 

• 2 31 

. 16* 

.4  66 

43 

.641  .2"4 

3 

. 89% 

1. 324 

1.7  86 

40 

• 141  1.303 

4 

1. 825 

3. 2*0 

3.67? 

39 

.067  3.3  17 

61.44  3?. 

14  K26A27A 

1.600  .310 

.016 

2.130  .2 70 

.30  7 

? • 94  0 .76  0 

.614 

4.  1401.340 

1.013 

6.2801 . 736 

. 666 

w?6A?7AS 


V * 

6 SP°E  * 

32.  1 

4 

apparent 

DIP 

LAvcp  THICKNESS 

velocity 

AT  ORIGIN 

1 

1. 44 0 

0 

. 916 

2 

1.600 

3 

.3  07 

3 

?.  130 

0 

.616 

4 

2. <540 

0 

1.  013 

6 

4.  140 

0 

. 866 

t> 

6.  ?«0 

0 

H 

OEPTM  A 

rMTCK  A 

VELOCITY  THICK  3 

DEPTH  4 

1 

.016 

.016 

1.440 

9 16 

.01  6 

2 

.323 

.3  07 

1.630 

307 

.72  3 

3 

.439 

• 616 

? • 1 30 

6 16 

.939 

4 

1.962 

1.  313 

2.990  1. 

on 

1 • °9? 

6 

2.413 

• 866 

4.140 

866 

2.818 

6 

6.28 

LAVER  ray  tc 

LAYER 

CRITICAL  31 S T 

RAY 

PROP  LAYER 

X 

Y 

X Y 

0.00 

0.00 

32.  14  0.00 

1 

.033 

.016 

.0  66 

32 

.107  .016 

2 

. 3 64 

.32  3 

.728 

31 

.776  .323 

3 

.824 

.0  31 

1.667 

31 

.311  .439 

4 

1. 662 

1.46? 

3.124 

30 

.679  1.94? 

6 

2.164 

2.818 

4.329 

29 

.476  2.810 

ip- 


75 


ei. 

44  34.44 

*263? 74 

t.hOf  . 

010 

.016 

2.100  . 

360 

. 424 

3 . 4 3 C i . 

060 

. 767 

4.4311. 

406 

. 969 

L 

4.  AM?. 

930  3 

.243 

1- 

W26H2  74 S 

N * 

6 SPRfAO  ■ 34 

• 6 

N 

APPARENT 

nip  catfr  tmtcknfss 

VELOCITY 

Af  ORIGIN 

1 

1.44V 

0 

. 016 

i 

2 

1 .600 

a 

• 4?4 

3 

2.  100 

0 

. 76  7 

4 

J.  4 J 0 

0 

. 94  9 

l* 

S 

4.43V 

a 

3.  ?4  3 

r» 

6 

6.  6 AO 

0 

■ 

N 

OFPTM  a 

THICK  a 

YTLOCITY 

THICK  « 

OF  PTH  « 

1 

• 3 1* 

.014 

1.440 

. 0 16 

.016  r 

2 

• 440 

• 4 ?4 

1.600 

.4  ;4 

• 44  0 

3 

U 197 

. 747 

?•  too 

.767 

1.147 

4 

?.  166 

. 944 

1.43  0 

. 96  4 

2.  146 

6 

4.399 

3*  243 

4.430 

3.243 

6.349 

6 

6.09 

\ 

LAYfR  RAY  TO 

HYFR  CRITICAL  CIST  RAY 

nro r layfr 

X 

Y 

t Y 

9.  00 

0. 00 

34.44  1.00  » 

f 

. 0JJ 

.01* 

.0*6 

34 

.407  .016  f 

| 

.414 

.440 

1 .029 

33 

.426  .440  1 

S 

.017 

t .147 

1 .634 

33 

.623  1.147 

4 

t.  741 

?.  1*4 

1.9ft 

32 

• 6V  9 2 • l 4 h 

S 

3.628 

4. 344 

r ,?66 

30 

.912  6.349 

VI 

PROGRAM  TO  COMPUTE  CORRECTIONS  TO 
HELICOPTER-RECORDED  DATA  (BEAUHELI ) 

AND  INPUT  AND  OUTPUT  DATA  FOR  PROCESSED  LINES 


% 


77 


PROGRAM  °E  A UHE L I 
P 1=3, t VI 592653  6 
OTP=PI/l «0. 

RTD*1«0.  /PI 
READ(5,100IAID1,AID2 
1 1 0 F OR  MAT  ! 2 AS ) 

REA  0(5, 105) IYRHOl,IOAl*lHRlf MINI ,01  FI , ITRM02, 

1 I0A2* IHR?,HIN2,0IF? 

IFC IYRH01.ME.IYRH02TG0  TO  5 
105  FORMAT (2 ( I V , 12  * X f 21 2 , X , F6. 2 ) ) 

GO  TO  10 

5 NRITE(6t*1101IfRMOl,lYRM02 
113  FORPATC*  CLOCK  CORR  YM  X,I4,X  OIFFEQS  FROM  YM  X,T4) 

CALL  EXIT 

10  T IH1  = I DA  1*56  20  0 ♦ IHR1 • 36 00 ♦ HI Nl*60 
TIM?^rOA2*562OO#IMR2*J60O*MlN2*6O 
RE A C (5 , 1 15 ) IS9YM, ISBO,IS8H,IS8P,ISBN,rSBP 
115  FORMAT (14, 12  *X , 2 12  »X ,212) 

S9TI**IS°0**6200»IS8h*3600*TSB*#68 
RFAC(5,120)CLK,SLK, WATVELK,08®,09E,TDtST,OR,0?R 
120  FORMAT I8F8. 3) 

WRITE!*, 125) AIDl,AI02 
125  FOR* AT ( X < , 2 A8 ) 

WRITE (6, 130) IYRM01,I0A1*IHR1,MIN1,DTF1 
130  r ORMAT ( X PRESHOT  TIME  CHECK  X , 1 4 , 12, X ,2 12 , X DELTA  T=  X, 

1 F7.JI 

WRITE  (6,  135)  IYRH02,  It)A2*lHR2,MlM2,0IF2 
135  F ORMAT ( X POSTSHOT  TIME  CHECK  X , 14 , 1 2 , * , 21 2 • 

1 X DELTA  T»  X,F 7.31 

WRITEI6,  140) 

140  F ORMAT ( X ALL  DISTANCES  IN  KILOMETERS,  ALL  VELOCITIES  X, 

1 XI*  KN/SEC,  ALL  T I W€  S ARE  UNIVERSAL  T IME*) 

WRITE  (6,  1451CL<,SLK,  WATVELK 
145  FORMAT (X  CARLE  LENGTH*  *,F*.3,X  ACTIVE  SECTION  X, 
lXLENGTHs  X ,F8«  3 , X WATER  VELOCITY*  X,F«.3) 

WRITE (6, 14  7) 098 ,D8E ,T  01  ST, OR  ,OBR 
147  FORMAT! t START  DEPTH*  X, 

1 F8.jf*  ENO  OERTH*  X,F8.3,X  LINE  LENGTH*  x,F8.3, 

2 x RECEIVER  OERTM=  X,FS.3,X  DEPTH  AT  RECEIVE**  X,F8.3) 
READ!4,15  0I NYRHOl.NCAl ,Nh©i, NHN1,NLAT1  , AL A T *1 , NL ONI , AL C N*1 

15  8 F0RMAT(I4,I2,X,2I2,X,  I3,X,F5.2,X,I4,X,F«5.2) 
ANTIH1=NOA1*86200»NHR 1* 36 00 ♦NNNl* 6 0 
ALATl*(NLATl*ALATHl/60.)*DTR 
ALON1=(NLON1-ALOMM1/60.V*DTR 
<TT  = 0 

IS  REA  0(4, 150) NTRM02,NrA2,NHR2,NPN2 ,NL AT  2 ♦ AL AT  M2 f NL ON2 , AL CN*2 
ANTIM2*NOA2*86200*NHR2*3600^NHN2®60 
ALAT2*!NLAT2*AL AT*?/E8. )*DT© 

ALON2=(NLON2-ALONM2/60. ) *DT9 
▼F(KIT.EQ.I)  GO  TO  28 

IF !S0T IH.GE.ANTI M1.AN0,S8TI*.LT. ANTIM2  ) GO  TO  20 
A NT  I Ml = A NT  I M? 

ALAT1*AL  AT2 
ALON1* Ai.  ON 2 
GO  TO  15 

29  FACT0P*!S8TIH- ANTIMll /( ANTTM2-ANTI*1) 

S9L0N*(AL0N2-AL0N1) •FACTORS  A LON1 
SBL AT= (ALAT2-AL A T1)*F ACTORS AL ATI 
S8LOP*SBLON*RTD 
SBLAP*SBLAT*RTO 

WRITE! 6, 155) ISBN , ISBP  » I S3Y *, ISBO, I SPH , ISB *, SELAP  ,S?LOP 
155  FORMAT ( X SONOBUOY  X,I2,<  AT  POSITION  X, 12, 


t t TIHF  *, 14 ,I2.X ,2I29<  L AT  X.F7.3,*  LONG  *» FP.3) 

PR  I TF  (6*  160> 

161  FORMAT  ( t OTRNOOA  shot  latitude  longitude  size  ^U«N  r, 

I * DEPTH  SO  I ST  CAOIS  TP  A VT  TTRAr  SURF  flCTC  t . 

Z FCLTC  SCOR  SHOTTIHF  SP*FL  SPVELC*) 

KlT*t 

?S  R£AO*Sf  i 651  I SHOT  .S  I ZT  . BURN  . OF  PK  , C I S T , K SY  H,  « SO  A « 
i XSHR*KS“N, SSEC .KHHR. ¥HHN,HSEC  .SPVEL 
I F | £ OF  ( S ) ) CALL  EXIT 

166  FOR-AT  tI*,4FS.? ,t,I4. I? .X,?l?,X,FS.2tX  ,21?, X.FS.2,  xfF6.5> 
SHORT  I n*«S0A*S6  2 30  ♦XSHR  *3  60  3 ♦« 5MN*6 0» SSEC-R U* * 

28  IF{ SHORT  IN. GE.ANTIHI.ANO. SHOSTIP.LT .ANT IP21G0  TO  10 
ANT  I N1 = A NT  I N? 

ALATt*AL at? 

ALONIsAL  on? 

GO  to  !«; 

30  FACTOR-1 SHCRTI  N-ANT I - l)  /( A NT  IN?- ANT  IN  1 ) 

SML  ^n*<alon?-aloni >*f  actor ♦aloni 

SHL  AT^ | ALAT2-AL  AT  11 »F  ACTOR»AlA7 l 
ARGsSlNF (ALAT1I*SINFI ALAT21  > 

1 COSFt  AL  ATI)  *COSr(  ALAT2)*C0SF(  ALON2-AL  ON  1 ) 
pOIST»636R.96R*ACOSFCIRG) 

ARG*SINF  C SOL  AT  ) • S I NF  T < NLA  T > ♦ 

1 COSFISBLAT)  •COSF<SHLAT)*COSF(  SHL0N-S9LCN7 
SOISTsS3^R.OSq*ACOSF«  ARC) 

SPVFLC=POI ST/f  ANTI N?- ANT  INI) 

SHOTX-SP  VEL*PURN-SLK-CL  X 
SHOTZ*  (-4Q.SSM.  St  J*«URN>  / 1 3 0 0. 

I c C SHOT? .GT. OEPX) SHOT  Z*OEPK 

SHf.OR-»SORT  < SHOTt  •SHOTX»SHOTZ*SHOT2)  /nATVE\K 

ST  !NC*SSFC  -SHCOR 

SHT  I NC=SHORTIH»  HORN- SHCOR 

CKCOR*OI Ft ♦ (0IF2-0*F !>• tSMCRTIP-f I PI ) / ( T IN? -T I N1 ) 

ART IHH»xSDA*«62  00»KHH»*36 J8»  XHPN*F0 ♦HSEC*C»COR 

TT«>AW«ARTtNM-SHT  INC 

TDEPsOPP  ♦«  OBF-OPH)  •OIST/THI  ST 

SI*  SORT  ( l.-PATVELK*NATm*/S./S.| 

SJRFCO*P=  I DR  ♦SHOT  Z)  •Sl/NATVfLX 
8 0TC0R=f  T3£»-oEPX  *08? -CPR) •St/WATVELK 
TTCOR*  Ttoa  N*SU»FC0R4 “OTCOR 
SHL  ON*SHLON*OTO 
LLON*$HLON 

ALONHa-(  SHLCN-LL0N)*6«. 

SHL  AT*SmL AT*RT  0 
LL  A T*SHL  AT 

ALATNafSHLAT-LLAT)*60. 

pRITE(6,  tZniKSTHfKSOA , i shot  « ll a t , al  at  n.llon , al onn, 

1 SI  7E*PI'RN,PF»X,S0IST  ,OIST,  TTCCRt  TT RA R , SURFC OR .POTCOR. 

? C*COR,SHCO»,KSHR, KSNK, ST  INC  * SPVEL  t SP VFLC 

WRITF<7,  17S|KSYP,<S0A  ,ISH0T,LLAT  ,ALATN,LLON,ALCNMf 

1 S I ZE,  °l'RN,  OF°Xf  SO  1ST  ,0 1ST.  T TC  FR  . T TR  A W , S L'RFCP  * . POT  C PR, 

? 0KC0R.SHC0O.KSHR ,KSNNt  ST IHC . SPVEL • S°VrLC 
170  FORMAT ( t TVI4. I?.X,I%vx.I3fXvFS.2.X.IV.Y*FS.?9XvFS.l.F6«l« 
l F6.J9pF6.?.6FS.2»X#2I2  »F6#?»?Ff  .S) 
t?«  FORPATC  l4,I2.X#I4,X913tX.FS.  ?*X,  m.X.FS.?*  XfFS.l,FF.l  , 
i F6.3.4F6. ?.4FS.2,X,2I?,F6.?92f?.S> 

GO  TO  25 
F NO 


UA[A  tHPUl  fUKMnJ  KJK  PKOGKAM  HLAUHttl 


IWHk.KH^I.HSfcC  . Hour  ,#»loute  and  Second  of  the  helicopter  ground  arrival. 

S^VE  L : Ship  velocity. 
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NAVIGATION  INPUT  FORMAT  FOR  PROGRAM  BEAUNEU 
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Line  identification. 
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Position  in  longitude  ideqrees  and  minutes)  of  the  snip's  navigation  points. 
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Figure  8.  Line  GL-1  record  section  and  velocity-depth  model 
interpretation.  Below  each  seismogram  is  listed  the  shot  number 
and  the  charge  weight  in  pounds. 
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Figure  9. 


Line  GL-2  record  section  and  velocity-depth  model  interpretation. 
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Figure  13.  Line  4-5  record  section  and  velocity-depth  model  interpretation 
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Figure  14.  Line  8A-9A  record  section  and  velocity-depth  model 
interpretation. 
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Figure  16.  Line  10-11  record  section  and  velocity-depth  model 
interpretation. 
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Figure  22.  Line  24B-25  record  section  and  velocity-depth  model  interpretation. 
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Figure  24.  Line  22B-23  record  section  and  velocity-depth  model  interpretation. 
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Figure  25.  Line  26A-27A  record  section  and  velocity-depth  model 
interpretation. 


-A/ 5 (SECONDS)  DEPTH 


DISTANCE  ( KM  ) 


6.88  KM/S 


DISTANCE  (SECONDS) 


Figure  26.  Line  26B-27B  record  section  and  velocity-depth  model 
in  terpre  tation . 
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